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This  report  describes  the  finite  element  computer  program  SOILSTRUCT,  used  in 
the  evaluation  of  soil-to-structure  interaction  of  earth  retaining  structures.  The  initial 
version  of  the  program  was  developed  by  Professors  G  W.  Clough  and  J.  M.  Duncan  in 
1969  and  has  been  enhanced  during  the  last  20  years  by  P.  »or  Clough  ana  his  co¬ 
workers.  This  report  documents  the  version  of  the  program  that  reflects  the  modifica¬ 
tions  made  in  conjunction  with  project  work  performed  at  the  US  Army  Engineer 
Waterways  Experiment  Station  (WES).  This  work  has  been  funded  by  the  Numerical 
Model  Maintenance  Program. 
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of  Mr.  Paul  Senter,  Chief,  CAED,  and  Dr.  N.  Radhakrishnan,  Chief,  ITL.  This 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric)  units 


as  follows: 

Multiply 

By 

To  Obtain 

cubic  feet 

0.2831683 

cubic  metres 

feet 

0.3048 

metres 

pounds  (mass)  per 

16.01846 

kilograms  per 

cubic  foot 

cubic  metre 

pounds  (force)  per 

47.88026 

pascals 

square  foot 

pounds  (force)  per 

6.894757 

kilopascals 

square  inch 

square  feet 

0.09290304 

square  metres 

square  inches 

6.4316 

square  centime 

/ 
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USER’S  GUIDE  FOR  THE  INCREMENTAL 
CONSTRUCTION  SOIL-STRUCTURE  INTERACTION 
PROGRAM  SOILSTRUCT 

PART  I:  DESCRIPTION  OF  PROGRAM  SOILSTRUCT 

Introductiou 

1.  SOILSTRUCT  is  a  general-purpose,  finite  element  program  for  two-dimensic 
plane  strain  analysis  of  soil-structure  interaction  and  soil-inclusion  interaction 
problems.  It  calculates  displacements  and  stresses  due  to  incremental  construction 
and/or  load  application  and  is  capable  of  modeling  nonlinear  stress-strain  material  be¬ 
havior.  The  simulation  of  incremental  construction  may  include  embankment  construc¬ 
tion  or  backfilling,  the  placement  of  layer(s)  of  a  reinforcement  material  during 
backfilling  or  embankment  construction,  dewatering,  excavation,  installation  of  a  strut 
or  tie-back  anchor  excavation  support  system,  removal  of  the  same  system,  and  the 
placement  of  concrete  or  other  construction  materials.  The  incremental  loading  simula¬ 
tion  may  consist  of  the  application  of  concentrated  loads,  boundary  pressures,  or  loads 
due  to  temperature  changes  in  non-soil  materials. 

2.  The  initial  version  of  SOILSTRUCT  was  developed  by  Professors  G.  W.  Clough 
and  J.  M.  Duncan  for  use  in  the  analysis  of  Port  Allen  and  Old  River  U-frame  locks 
(Clough  and  Duncan  1969).  This  version  of  the  program  reflects  modifications  made  in 
conjunction  with  a  number  of  projects  at  the  US  Army  Engineer  Waterways  Experiment 
Station  to  expand  the  capabilities  of  the  finite  elements  constitutive  models,  load  vector 
formulation  algorithms,  the  size  of  the  problem  which  may  be  analyzed,  and  the  transfer 
of  input,  output,  restart,  and  plot  data  files  by  means  of  disc  storage.  SOILSTRUCT  has 
been  coded  in  Fortran  77  language  and  consists  of  a  main  program  and  26  subroutines 
named:  DETNA,  INITAL,  STRSTF,  QUAD,  BAREL,  EXCAV,  EQNDFO,  3URFLD, 
JTSTF,  SUBSTP,  JSTRES,  SEEP,  MODCAL,  BUILD,  OPTSOL,  AUXOUT,  STRESS, 
PRNCIP,  PRNTFD,  GETFIL,  CNVERT,  REBAR,  FILLBARS,  BARSTIF,  RESTRESS, 
and  NFACTS.  A  user’s  guide  for  the  SOILSTRUCT  Program  can  be  found  in  Appen¬ 
dix  A.  Appendix  B  contains  the  sequence  of  operations  for  the  SOILSTRUCT  Program. 
Appendix  C  applies  simple  beam  theory  and  the  theory  of  elasticity  to  the  bending  of 
structural  members. 


Finite  Elements  Employed 

3.  Three  types  of  finite  elements  are  used  to  represent  the  behavior  of  different 
materials:  (a)  a  two-dimensional  continua  element,  (b)  an  interface  element,  and  (c)  a 
one-dimensional  bar  element. 

4.  A  two-dimensional,  subparametric,  quadrilateral  element  (QM5)  is  used  to  repre¬ 
sent  the  soil  and  most  structural  materials.  Structural  supports,  such  as  the  struts  or 
tieback  components  of  an  excavation  support  system,  are  typically  modeled  as  a  spring 
support  using  bar  elements.  However,  two-dimensional  elements  have  been  used  to 
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model  these  supports.  The  geometry  of  this  element,  developed  by  Doherty,  Wilson, 
and  Taylor  (1969),  is  defined  by  four  external  nodes,  while  the  displacement  functions 
include  an  internal  fifth  node.  To  improve  flexural  response,  a  constant  shear  strain, 
calculated  at  the  location  of  the  interna',  fifth  node,  is  imposed  throughout  the  element. 
The  QM5  element  can  be  allowed  to  degrade  to  a  triangular  element  by  letting  two  ad¬ 
jacent  nodes  of  the  quadrilateral  coincide. 

5.  The  Goodman,  Taylor,  and  Breeke  (1968)  interface  element  is  used  to  allow  for 
relative  movement  between  different  materials,  such  as  between  a  soil  backfill  and  a 
support  wall.  This  element  is  defined  by  four  nodesi,  with  each  of  the  two  pairs  of 
nodes  having  the  same  coordinates;  thus,  this  type  of  element  has  no  thickness. 

6.  One-dimensional,  two-node,  bar  or  spring  elements  are  used  to  model  the  be¬ 
havior  of  a  variety  of  structural  systems.  This  includes  the  modeling  of  structural  sup¬ 
ports  such  as  braces  or  tiebacks  or  the  modeling  of  reinforcement  placed  within  a  soil 
backfill. 


PART  II:  MATERIAL  STRESS-STRAIN  BEHAVIOR 

7.  Several  inodes  of  stress-strain  behavior  are  utilized  to  represent  the  response  of 
soil,  construction  materials,  and  the  inteiface  region  between  different  materials. 

8.  The  constitutive  relationship  used  for  all  two-dimensional  elements  is  Hooke’s 
law.  SOILSTRUCT  uses  an  incremental,  equivalent  linear  method  of  analysis  to  model 
nonlinear  material  behavior.  In  this  type  of  analysis,  the  incremental  changes  in  stress¬ 
es  are  related  to  the  incremental  strains  through  a  linear  relationship.  This  relationship 
is  defined  for  each  structural  element  by  two  engineering  constants,  the  Young’s 
moduli  and  the  Poisson’s  ratio.  For  the  soil  elements,  either  the  Young’s  moduli  and 
Poisson’s  ratio,  or  the  Young’s  moduli  and  bulk  moduli  may  be  specified. 

Nonlinear  Stress-Strain  Response  of  Soil 

9.  A  plane  strain,  isotropic  drained  or  undrained  stress-strain  soil  model  is  incor¬ 
porated  within  SOILSTRUCT.  The  program  uses  a  nonlinear,  stiess-dependent  hyper¬ 
bolic  curve  to  represent  the  relationship  between  stress  and  strains  developing  during 
primary  loading  of  the  soil  (Figure  la)  and  a  linear  stress-strain  response  during  unload¬ 
ing  or  reloading  of  the  soil  (Figure  lb).  The  unload-reload  stress-strain  response  is 
applicable  when  the  current  stress  state  is  less  than  that  which  has  been  applied  pre¬ 
viously;  otherwise,  the  primary'  loading  stress-strain  is  appropriate.  Laboratory  testing 
and  interpretation  procedures  for  determining  the  parameters  used  to  define  the  soil 
model  are  described  in  Duncan,  Byrne,  Wong,  and  Mabry  (1978).  A  brief  review  of  the 
hyperbolic  model  is  given  in  paragraphs  12-17. 

10.  The  nonlinear  soil  response  to  loading  is  modeled  by  performing  a  series  of 
analyses  in  which  each  load  is  applied  incrementally,  with  the  total  change  in  stress 
computed  at  the  center  of  each  soil  element  being  equal  to  the  sum  of  the  incremental 
changes  in  stress  over  all  the  load  steps.  In  general,  the  greater  the  curvature  of  the 
stress-strain  relationship  or  the  larger  the  magnitude  of  the  applied  load,  the  greater 
the  number  of  load  steps  required  to  accurately  model  the  nonlinear  soil  response. 

This  may  be  achieved  in  two  ways  using  SOILSTRUCT;  either  the  total  load  is  ap¬ 
plied  using  a  greater  number  of  incremental  loadings,  or  during  the  course  of  each  load 
case  analysis,  the  load  vector  may  be  applied  in  a  series  of  increments  using  the  sub¬ 
step  option. 

1 1.  Application  of  each  loading  in  the  finite  element  analysis  results  in  a  change  in 
stress  within  each  of  the  soil  elements.  In  addition  to  the  change  in  stress,  there  is  a  cor¬ 
responding  change  in  stiffness.  Since  each  incremental  analysis  is  performed  assuming 
equivalent  linear  element  response,  SOILSTRUCT  updates  the  value  of  the  elastic 
moduli  assigned  to  each  soil  element  so  as  to  reflect  the  magnitude  of  the  current  stress 
state  within  the  element.  To  account  for  the  change  in  stiffness  that  occurs  during  the 
application  of  a  load  increment,  each  incremental  load  calculation  may  be  repeated 
using  the  iteration  option.  When  the  iteration  option  is  invoked,  the  load  vector  is  reap¬ 
plied  with  a  revised  value  for  the  element  stiffness.  The  value  assigned  for  the  stiffness 
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a.  Hyperbolic  representation  of  stress-strain  curve  for  primary  loading 


b.  Linear  unloading-reloading  stress-strain  relationship 

Figure  1.  Hyperbolic  model  for  stress-strain  behavior 
(after  Duncan,  Byrne,  Wong,  and  Mabry  1978) 
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of  the  soil  element  reflects  the  average  of  the  stiess  state  developing  at  the  end  of  the 
previous  load  case,  or  substep,  and  that  which  develops  during  the  current  iteration. 
However,  when  only  one  iteration  is  specified,  the  modulus  values  am  calculated  usin^ 
the  stresses  developing  at  the  end  of  the  previous  load  increment.  Upon  completion  of 
the  last  iteration  for  each  load  case  or  substep,  the  arravs  tabulating-the  values  of  the 
total  nod’l  point  displacements  and  total  element  su —  updated  with  the  com¬ 
puted  incremental  values. 

Primary  Loading  -  Young’s  Moduli 

12.  Prior  to  each  analysis  a  tangent  Young’s  modulus  Et*  is  assigned  to  each  .1 
element.  The  stress-dependent  value  of  E(  is  computed  using  the  relationship 


where 

Ej  =  initial  Young’s  modulus 
Rf  =  failure  ratio 
SL  ■  s- /css  level 

The  initial  Young’s  modulus  Ej  is  equal  to 


where 

K  a  modulus  number 
Pa  =  atmospheric  pressure 
n  =  modulus  exponent 
O3  »  minor  principal  stress 

13.  The  proportion  of  mobilized  shear  strength  for  each  soil  element  is  reflected  in 
the  value  of  the  -tress  level,  SL.  SL  is  equal  to  the  current  deviator  stress  (Oj  -  O3) 

divided  by  the  deviator  stress  at  failure  (Oj  -  o3)f  ,  denoted  by  the  subscript  f . 


For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation  (Appendix  D). 
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SL  = 


°l-q3 


(°1 "  °3)f  (3) 

where  o j  =  major  principal  stress.  The  value  of  SL  ranges  from  a  value  equal  to  zero 

to  a  value  equal  to  unity.  SL  equal  to  zero  indicates  an  isotropic  stress  state,  while  SL 
equal  to  unity  corresponds  to  the  complete  mobilization  of  shear  resistance  within  the 
soil  element. 

14.  This  version  of  SOILSTRUCT  allows  for  two  procedures  for  defining  the 
deviator  stress  at  failure;  the  original  Duncan  formulation,  Duncan  and  Chang  (1970), 
and  a  procedure  developed  by  Peters.  In  the  original  Duncan  formulation  as  shown  in 
Figure  2a,  the  value  of  the  minor  principal  stress  at  failure  is  set  equal  to  the  current 
minor  principal  stress.  The  deviator  stress  at  failure  is  given  by 

2c  cos  tp  +  2o,  sin  q> 

1  3  f  1  -  sin  9  (4) 

where 

C  =  cohesion  intercept 
9  »  angle  of  internal  friction 

In  the  Peters  formulation  shown  in  Figure  2b,  the  value  for  the  deviator  stress  at  failure 
is  defined  based  upon  the  assumption  that  the  average  value  of  the  major  and  minor 
principal  stresses  at  failure  is  equal  to  the  average  value  of  the  current  major  and  minor 
principal  stresses. 


(Of i  -  o3)f  =  2c  cos  9  +  (Cj+  o3)  sin  tp 


15.  The  difference  between  these  two  procedures  may  be  expressed  in  terms  of  the 
resulting  vector  curves  as  shown  in  Figure  2.  Vector  curves  are  loci  of  points  describ¬ 
ing  the  state  of  stress  on  planes  on  which  failure  eventually  occurs.  The  resulting  vec¬ 
tor  curves  in  the  Duncan  formulation  are  to  the  right,  reflecting  the  assumption  that 
failure  is  a  result  of  an  increase  in  major  principal  stress.  In  contrast,  the  vector  curves 
to  the  left  in  the  Peters  formulation  are  attributed  to  a  coincident  increase  in  major  prin¬ 
cipal  stress  and  a  decrease  in  the  minor  principal  stress.  In  general,  the  Duncan  for¬ 
mulation  results  in  larger  values  of  (Oj  -  o3)f  and  therefore  smaller  values  of  SL 

than  the  Peters  formulation.  The  Peters  formulation  was  developed  for  undrained 
loading  where  (Oj  +  o3)  is  the  constant  consolidation  stress  determined  for  the 
initial  stress  state. 


a.  J uncan  formulation,  a.  =  a 

3f  : 


•  0  »  .  .  <0]  +  o3)f  (Oj  +  03) 

b.  Peters  formulation,  p  -  — ! - *-±-  »  — i - i- 

2  2 


Figure  2,  Two  procedures  used  to  define  deviator  stress 
at  failure,  (a.  -  a,). 


16.  The  failure  ratio  Rf  relates  the  ultimate  deviator  stress  (Oj  -  to  the 

deviator  stress  at  failure  (Oj  -  Oj)f  . 


(6) 


The  ultimate  deviator  stress  is  the  asymptote  to  the  stress-strain  hyperbola,  as  shown  in 
Figure  la.  The  value  of  Rf  is  always  less  than  unity  and  varies  from  O.S  to  0.9  for 
most  soils. 


Unload-Reload  Stress-Strain  Behavior  •  Young’s  Modu'us 

17.  During  unloading  or  reloading,  when  the  current  deviator  stress  is  less  than  that 
which  has  been  applied  during  previous  loadings,  a  stress-dependent,  linear  response  is 
assumed,  as  shown  in  Figure  lb.  In  this  case,  the  value  of  Eur  is  computed  using 


Va 


(7) 


where  Kur  =  the  unload-reload  modulus  number. 

Poisson’s  Ratio  or  Bulk  Modulus 

18.  The  second  elastic  parameter  used  to  define  the  material  behavior  of  soil  is 
either  the  Poisson’s  ratio,  v  ,  or  the  bulk  modulus,  B  .  This  version  of  SOILSTRUCT 
allows  either  parameter  to  be  used.  When  using  Poisson’s  ratio,  two  values  are  speci¬ 
fied:  a  constant  value  which  is  applicable  for  all  states  of  stress  prior  to  failure,  SL  <  1 
,  and  the  value  of  Poisson’s  ratio  applicable  when  the  shear  strength  of  the  soil  is  fully 
mobilized,  SL  =  1  . 

Poisson’s  ratio  formulation 

19.  The  hyperbolic  model  is  designed  so  that  as  the  soil  approaches  failure,  v  =  0.5  . 
The  variation  in  v  is  accomplished  by  computing  the  shear  modulus  with  Et  and  Vj 
and  bulk  modulus  with  Ej  and  Vj ,  where  Vj  is  the  value  of  Poisson’s  ratio  before 
failure.  This  variation  in  v  amounts  to 


(8) 
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where 


iz  h 

1  +  V: 


Therefore,  v  *  v-  when  Et  =  Ej ,  but  increases  toward  O.S  as  Et  becomes  small  near 
failure. 

Bulk  modulus  formulation 

20.  Because  many  soils  exhibit  nonlinear  and  stress-dependent  volume  change  char¬ 
acteristics,  the  stress-dependent  bulk  modulus  formulation  developed  by  Duncan, 
Byrne,  Wong,  and  Mabry  (1978)  is  also  included.  According  to  the  theory  of  elasticity 
the  value  of  the  bulk  modulus  B  is  defined  as  the  ratio  of  the  change  in  mean  principal 
stress  to  the  change  in  the  volumetric  strain. 


AOj  +  Ao2  +  Ao3 


where  Actj  ,  A o2  ,  and  A03  are  the  changes  in  the  values  of  principal  stress  and  is 
the  corresponding  change  in  volumetric  strain. 

21.  The  bulk  modulus  for  soil  has  been  found  to  increase  in  value  with  increasing 
values  of  minor  principal  stress  and  is  assumed  to  be  independent  of  stress  level.  It  is 
approximated  by  the  equation 


B-  KbPa 


where 


bulk  modulus  number 


m  =*  bulk  modulus  exponent 

22.  Experience  with  this  formulation,  as  described  by  Duncan,  Seed,  Wong,  and 
Ozawa  (1984),  has  led  to  the  following  restrictions  on  the  value  assigned  to  the  bulk 
moduli.  The  maximum  value  of  bulk  moduli  corresponds  to  the  situation  when  the 
value  for  Poisson’s  ratio  v  is  equal  to  0.49.  The  value  of  bulk  moduli  may  be  defined 
by  the  values  of  Young’s  modulus  Et  and  v  as  given  by  the  relationship 


B-  A  I— I— 1 
3  L(l  -  2v)  J 
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(12) 


The  bulk  moduli  attains  a  maximum  value  when  v  =  0.49  and  is  equal  to 

(13) 

23.  The  value  assigned  for  the  bulk  modulus  is  also  restricted  to  a  minimum  value. 
This  minimum  value  is  determined  from  the  relationships  between  bulk  modulus, 
Poisson’s  ratio,  the  at-rest  earth  pressure  coefficient  K0  ,  and  the  angle  of  internal  fric¬ 
tion  <p .  For  a  single  linear  elastic  material  (Dunlop,  Duncan,  and  Seed  1968)  KQ  and 
v  are  related  as  follows: 


v 


1  +  K0 


(14) 


Jaky  (1948)  suggested  that  K  may  be  approximated  using  the  relationship 


Kq  =  1  -  sin<p 


By  introducing  Equations  14  and  15  into  Equation  12,  the  minimum  value  for  the  bulk 
moduli  value  during  primary  loading  is  restricted  to 


Experience  gained  when  using  this  soil  model  to  simulate  excavation  resulted  in  the 
development  of  the  restriction  on  the  value  of  bulk  moduli  during  unloading-reloading 
behavior  as  given  by 


This  relationship  was  obtained  by  introducing  Equation  14  into  Equation  12. 

Consolidation  stress  method  for  undrained  analysis 

24.  Soil  strength  is  typically  defined  in  SOILSTRUCT  by  cohesion  c  and  friction 
angle  9 ,  which  are  chosen  to  be  appropriate  for  the  drainage  condition  of  each  ele¬ 
ment  based  on  its  permeability  and  the  loading  rate.  For  undrained  conditions,  how- 
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ever,  this  approach  is  not  suitable.  In  order  to  model  the  increase  in  strength  produced 
by  higher  consolidation  stress  it  is  necessary  to  either  assign  a  different  cohesion  (with 
<p  =  0)  to  each  element,  which  is  not  practical,  or  to  assign  a  total  stress  friction  angle 
to  each  material,  which  is  physically  inconsistent  for  saturated  materials.  The  correct 
result  can  only  be  obtained  for  undrained  conditions  by  selecting  the  undrained  strength 
from  the  preloading  consolidation  conditions  and  setting  cp  =  0  for  all  subsequent  un¬ 
drained  loadings.  Therefore,  the  program  allows  the  strength  to  be  input  as  a  ;-iio .... 
strength  to  effective  consolidation  pressure  (Su/pg  ) .  ise,  n  u  set  to  zero  and 

K  is  expressed  as  a  function  of  the  initial  consolidation  stress;  this  is  an  approach 
similar  to  that  described  below  for  soil  strength  whereby  the  strength  is  based  on  the  ini¬ 
tial  consolidation  state  and  the  friction  angle  is  set  to  zero.  It  has  been  found  through 
experience  that  the  initial  modulus,  Pa  K ,  can  be  expressed  as  a  ratio  of  the  undrained 
shear  strength  (Clough  and  Tsui  1977;  Mana  1978)  whereby  E;  =  KSU  .  Thus,  assum¬ 
ing  the  constant  K  is  known,  the  undrained  shear  strength  becomes  the  fundamental 
parameter  controlling  the  response  of  the  soil.  The  procedure  consists  of  the  following; 

a.  The  consolidation  stress  is  computed  for  each  element  based  on  the  geometry 
and  boundary  conditions  prior  to  loading  assuming  that  the  pre-existing  ele¬ 
ments  have  fully  consolidated  under  their  own  weight.  Elements  above  the 
water  table  are  assigned  the  total  unit  weight  of  the  soil,  and  elements  below 
the  water  table  are  assigned  the  buoyant  unit  weight.  The  stresses  created  by 
this  configuration  are  computed  by  subroutine  INITAL. 

b.  The  effective  consolidation  stress  pc  is  computed  for  each  element  as 


where  o'  and  o'  are,  respectively,  the  horizontal  and  vertical  effective 

stresses.  This  value  is  stored  for  each  element  for  use  in  all  subsequent  cal¬ 
culations. 

c.  Each  material  type  is  assigned  a  value  of  Su/p'  and  K  .  These  values  are 

then  combined  with  p'  computed  from  the  initial  stress  computations  to 

determine  Su  and  E;  for  each  element.  The  property  values  assigned  to 
each  element  therefore  depend  on  material  type  and  section  geometry.  For 
example,  shear  strengths  may  be  moderately  higher  under  a  levee  centerline 
than  at  the  toe  as  a  result  of  the  higher  consolidation  stress  imposed  by  the 
levee. 

Poisson’s  ratio  for  initial  stress  computation 

25.  When  computing  initial  stresses  by  gravity  turn-on  analysis,  the  value  of 
Poisson’s  ratio  used  for  the  soil  model  may  not  be  suitable  for  initial  stress  conditions. 
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For  example,  the  Poisson’s  ratio  for  undrained  analysis  is  generally  taken  to  be  nearly 
0.5  because  saturated  soil  is  nearly  incompressible.  By  contrast,  the  initial  stress  condi¬ 
tions  should  be  fully  drained.  Therefore,  the  initial  stress  computations  are  based  on  a 
value  of  Poisson’s  ratio  that  gives  the  correct  ratio  o'  /o'  for  level  ground  condi¬ 
tions;  that  is,  o'  /o^  =  K0  .  Using  the  value  of  K0  input,  the  Poisson’s  ratio  used 
forinitial  conditions  is  given  by  Equation  14. 
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PART  IU:  MODELING  STRUCTURAL  ELEMENTS 
Structural  Material  Response 

26.  Structural  materials  such  as  wood,  concrete,  or  steel  are  modeled  using  two- 
dimensional  QM5  elements  with  linear  elastic  stress-strain  t^avior  assumed.  Support 
elements  such  as  struts  or  anchors  are  typically  modeled  using  one-dimensional  bar  ele¬ 
ments,  and  are  also  assumed  to  behave  linearly.  Bar  elements  as  formulated  within 
SOILSTRUCT  have  the  capability  to  respond  in  compression  only,  in  tension  only,  or 
in  both  tension  and  compression.  In  addition,  slack  in  the  support  system  at  the  time  of 
installation  may  be  accounted  for  by  specifying  an  initial  value  of  displacement  for  the 
bar  element. 

Soil  Reinforcement 

27.  Reinforcement  placed  within  a  soil  backfill  provides  resistance  to  tensile 
strains  which  may  develop  as  a  result  of  loadings  attributed  to  the  force  due  to 
gravity  acting  on  the  soil  mass  or  due  to  applied  loadings.  The  reinforcement  within 
a  soil  backfill  or  an  embankment  is  modeled  using  one-dimensional  bar  elements. 
Reinforcement  has  the  effect  of  increasing  the  stiffness  of  the  soil  mass  and  is 
modeled  by  increasing  the  stiffnesses  of  each  of  the  two-dimensional  soil  elements 
that  includes  a  layer  of  reinforcement. 

28.  Reinforcement  can  be  modeled  by  either  of  two  methods.  The  first  method  con¬ 
sists  of  using  bar  elements,  as  tension-only  elements,  sandwiched  between  solid  ele¬ 
ments.  By  this  method  the  location  of  reinforcing  must  be  accounted  for  in  laying  out 
the  analysis  mesh.  This  may  become  inconvenient  if  reinforcing  layers  are  closely 
spaced. 

29.  The  second  method  allows  the  reinforcing  to  be  “embedded”  into  the  interior  of 
the  element  such  that  it  is  unnecessary  to  place  reinforcement  at  element  boundaries. 

The  stiffness  of  the  reinforcement  is  added  to  the  solid  element  as  follows:  the  stiffness 
for  the  reinforcement  “bar”  element,  [Kb] ,  which  defines  the  relationship  between  the 
displacement  {ub}  and  the  bar  forces,  is  then  related  to  the  node  displacements  {u) 
by  the  interpolation  {ub}  ■  [N]{u}  ,  where  [N]  is  the  interpolation  matrix  for  displace¬ 
ment.  Thus  the  effect  of  the  bar  on  the  element  stiffness  is  obtained  by  adding  to  it  the 
quantity  [N]T[Kb][N].  The  procedure  can  be  repeated  for  any  number  of  bars  placed 

at  any  orientation  within  the  element.  Also,  the  procedure  does  not  add  to  the  total  num¬ 
ber  of  degrees  of  freedom.  The  only  input  required  to  add  reinforcement  by  method 
two  is  the  stiffness  and  location  of  the  reinforcement  layer;  the  individual  soil  elements 
affected  by  the  layer  are  determined  by  SOILSTRUCT. 

30.  During  each  analysis  the  reinforcement  model  monitors,  the  resulting  total  strains 
within  each  bar  element  and  distinguishes  between  tensile  and  compressive  strains. 
Typical  soil  reinforcement  materials  are  very  thin  relative  to  their  length  and  are  not 
capable  of  resisting  compressive  strains.  This  restriction  is  incorporated  within  the  rein¬ 
forcement  model  by  reducing  the  reinforcement  bar  stiffness  to  a  near  zero  value  when 
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compressive  strains  are  computed.  The  model  has  a  provision  for  re-establishing  the 
tensile  stiffness  of  the  reinforcement  bar  element(s)  to  the  two-dimensional  soil  ele¬ 
ments)  when  and  if  tensile  strains  occur  within  the  bar(s)  during  subsequent  loadings. 

Interface  Response 

31.  Interface  elements  are  used  to  allow  for  relative  movement  between  different 
material  regions,  such  as  between  a  soil  backfill  and  a  support  wall.  These  elements  are 
defined  by  four  nodes,  each  node  having  two  degrees  of  freedom;  each  of  the  two  pairs 
of  nodes  sharing  the  same  coordinates.  The  interface  element,  therefore,  is  of  finite 
length  but  zero  thickness. 

32.  The  properties  of  interface  elements  are  defined  by  an  interface  normal  stiffness 
kj,  and  an  interface  shear  stiffness,  ks  .  These  values  of  stiffness  relate  the  average 
relative  displacements  normal  to  the  interface  element  An  and  average  relative  shear 

displacements  A$  to  the  corresponding  normal  stress  on  and  shear  stress  t  by  the 
equations 


k  A 
n  n 


(19) 


and 


x  = 


k.  A 


s  s 


(20) 


The  units  of  kn  and  k$  are  force  per  cubic  length. 

33.  The  initial  value  of  kQ  is  set  equal  to  1  x  10s  within  the  program.  This  value 
for  kR  ensures  that  the  normal  relative  displacement  of  the  interface  element  is  insig¬ 
nificant  when  English  units  (*'et,  pounds)  or  SI  units  (meters,  kilonewtons)  are  used.  If 
other  units  are  used,  the  value  of  the  normal  stiffness  may  need  to  be  changed  to  a 
higher  value. 

34.  Two  types  of  interface  shear  response  are  modeled,  a  bilinear  shear  stress-dis¬ 

placement  relationship  shown  in  Figure  3a,  and  a  hyperbolic  shear  stress-displacement 
relationship  shown  in  Figure  3b.  In  the  bilinear  model,  the  value  assigned  to  ks  is  a 
constant  so  long  as  the  average  shear  stress  x  along  the  interface  is  less  than  the  shear 
strength.  If  the  shear  strength  of  the  interface  element  is  fully  mobilized,  which  occurs 
when  t  is  equal  to  x{ ,  k#  is  set  equal  to  zero.  When  the  normal  stress  is 

greater  than  or  equal  to  zero,  the  value  of  Xf  is  given  by  the  relationship 


RELATIVE  SHEAR  DISPLACEMENT,  As 


a.  Bilinear  stress-strain  model  representing  interface  behavior 


b.  Hyperbolic  representation  of  the  variation  of  shear  stress  with  relative 
shear  displacement  (kst  =  tangent  interface  shear  stiffness, 
ktf  »  initial  interface  shear  stiffness) 

Figure  3.  Bilinear  and  hyperbolic  models  for  interface  shear  stress-relative  shear 
displacement  behavior  after  Clough  and  Duncan  (1969) 
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xt  =  c.  +  a  tan  8 

tin 


(21) 


where  ' 

Cj  =  cohesion  intercept  along  the  interface 
5  =  angle  of  internal  friction  along  the  interface 
and  shown  in  Figure  4a.  When  on  is  less  than  zero,  x{  is  computed  using 


where  ot  =  tensile  strength. 

35.  Direct  shear  test  results  on  soil-to-concrete  interfaces  and  soil-to-steel  interfaces 
by  Potyondy  (1961),  Clough  and  Duncan  (1969),  and  Peterson  et  al.  (1976)  have  shown 
that  the  value  of  5  is  proportional  to  the  angle  of  internal  friction  of  the  soil.  The 
value  of  the  constant  of  proportionality  is  dependent  upon  both  the  type  of  soil  and  the 
type  of  material  comprising  the  surface  of  the  structure. 

36.  The  direct  shear  tests  performed  by  Clough  and  Duncan  (1969)  and  Peterson  et 
al.  (1976)  have  shown  that  for  some  materials,  such  as  sand-to-concrete  interfaces,  the 
interface  response  during  shear  is  nonlinear  and  dependent  upon  the  normal  stress.  A 
nonlinear,  stress-dependent  hyperbolic  curve  is  used  to  represent  the  relationship  be¬ 
tween  shear  stress  and  average  relative  shear  displacement  developing  during  primary 
loading  of  the  interface  (Figure  3a)  and  a  linear  shear  stress-relative  displacement 
response  during  unloading  or  reloading  of  the  interface.  The  stress-dependent  value  of 
kst  is  computed  using  the  relationship 


k,i<* 


(23) 


where 

ksj  a  initial  interface  shear  stiffness 
Rfi  a  failure  ratio 
SL{  a  stress  level 

The  initial  interface  shear  stiffness  ksi  is  equal  to 


(24) 


b.  Normal  stress-normal  relative  displacement  behavior 

Figure  4.  Interface  strength  criteria  and  normal  stress-normal 
relative  displacement  relationship 


where 


ksi  -  Vw 


\  J 


n. 


Kj  =  interface  modulus  number 
yw  =  unit  weight  of  water 
n;  =  interface  modulus  exponent 

37.  The  proportion  of  mobilized  shear  strength  for  each  interface  element  is 
reflected  in  the  value  of  the  stress  level  SL; .  SL.  is  equal  to  the  current  shear  stress 
x  divided  by  the  stress  at  failure,  xf . 


(25) 


Tf  is  computed  using  either  Equation  21  or  22.  SLj  ranges  in  value  between  zero  and 
one. 

38.  The  failure  ratio  Rf|  relates  the  ultimate  shear  stress  tult  to  the  shear  stress  at 
failure. 


Tf  *  Rfi  ’  Tult 


(26) 


The  ultimate  shear  stress  is  the  asymptote  to  the  shear  stress-relative  shear  displace¬ 
ment  hyperbola,  as  shown  in  Figure  4b.  Direct  shear  tests  on  sand-to-concrete  inter¬ 
faces  by  Peterson  et.  al  (1976)  have  shown  the  value  of  Rfi  typically  ranges  in  valve 
from  0.3  to  1.0. 

39.  The  relationship  between  the  average  normal  stress  along  the  interface  and  the 
tensile  strength  is  shown  in  Figure  4b.  The  value  of  kn  is  a  constant  value  equal  to 
l  x  108  when  on  is  greater  than  or  equal  to  ot.  If  on  is  Idss  than  ot,  kn  is  set 
equal  to  zero,  assuring  that  additional  tensile  stresses  do  not  accrue  upon  subsequent 
loadings.  This  procedure  allows  for  separation  to  occur  between  two  adjacent  regions 
of  the  mesh  along  interface  elements  during  the  course  of  an  incremental  analysis.  If 
the  separation  closes  during  subsequent  loading,  k„  is  reset  to  1  x  10®. 
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Sheet-Pile  Element 


40.  Due  to  the  slenderness  of  aheet  piles,  their  primary  mode  of  response  to  earth 
and  water  loadings  is  the  development  of  bending  stresses  within  the  sheet  pile.  Thus, 
the  response  of  the  finite  element  used  to  model  the  sheet  pile  in  bending  is  of  primary 
importance.  Representation  of  bending  stiffness  in  soil-strUw.-.*“  ^t^raction  analyses 
has  always  presented  a  difficulty.  If  an  element  is  formulated  for  bending  using  the 
approach  found  in  most  structural  analysis  codes,  an  incompatibility  is  created  between 
the  bending  and  solid  (soil)  elements.  This  incompatibility  results  from  the  technical 
requirement  that  displacement  gradients  (slope)  must  be  continuous  across  beam  ele¬ 
ments  whereas  the  solid  elements  generally  only  provide  for  continuous  displacements. 
The  incompatibility  problem  is  avoided  in  SOILSTRUCY  by  using  slender  solid  ele¬ 
ments  to  model  bending.  These  elements  are  similar  to  the  soil  elements,  rather  than 
true  beam  elements.  In  fact,  the  particular  choice  of  element  formulation  selected  for 
the  SOILSTRUCT  code  was  made  to  ensure  that  the  solid  elements  would  correctly 
model  strain  patterns  associated  with  bending.  Experience  by  Mana  (1978)  on  a  num¬ 
ber  of  soil-structure  interaction  problems  has  shown  this  approach  to  work  well. 

Sheet-Pile  Section  Properties 

41.  The  properties  of  the  solid  elements  used  to  model  the  sheet  pile  are  the  elastic 
properties,  E  and  v.  However,  the  stiffness  of  flexural  members  depends  on  the 
product  El  where  I  is  the  moment  of  inertia,  a  geometric  property  of  the  member’s 
cross  section.  The  finite  element  has  a  rectangular  cross  section  having  a  unit  width  (in 
the  out-of-plane  direction)  and  a  height  that  depends  upon  how  the  analysis  mesh  is 
drawn.  For  example,  a  sheet-pile  wall  has  i  complex  cross-sectional  shape,  but  is  repre¬ 
sented  in  the  finite  element  mesh  as  a  unit-wide  rectangle  with  a  height  that  depends  on 
such  considerations  as  maintaining  an  aspect  ratio  of  the  finite  element  that  is  favorable 
from  a  standpoint  of  numerical  accuracy.  To  achieve  the  response  to  bending  that  is 
consistent  with  simple  beam  theory  (see  section  titled  "Simple  Beam  Theory”  in  Ap- 
pendic  C),  the  modulus  of  the  element  must  be  chosen  to  obtain  the  equivalent  flexural 
stiffness  as  specified  by  the  product  El ,  where  I  is  the  moment  of  inertia  per  unit 
width  of  sheet  pile  wall.  Therefore,  the  properties  of  the  sheet-pile  elements  are  deter¬ 
mined  such  that  the  section  stiffness  of  the  plane  strain  element  Hgl^l  -  v|)  (see 

section  titled  “Pure  Bending  of  Members  in  Plane  Strain”  in  Appendix  C)  matches  the 
El  of  the  sheet  pile.  Therefore,  the  sheet  pile  elements  obtain  proper  bending  stiffness 
when  assigned  the  modulus  given  by: 

Ee  -  f!(.  -  *  <27> 

c 

where 

Ee  =  equivalent  Young’s  modulus  for  pile  elements 
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E  «*  Young’s  modulus  for  the  sheet  pile 
i  =  moment  of  inertia  per  unit  width  of  sheet-pile  wall 
Ie  =  equivalent  moment  of  inertia  of  the  element  used  to  model  the  sheet  pile 


4  =  j  2 


(28) 


and 

be  =  unit  width  of  the  sheet-pile  wall  mesh 
he  *  width  of  the  element  used  to  model  the  sheet  pile 
ve  =  Poisson’s  ratio  of  the  element  used  to  model  the  sheet  pile 

The  factor  (1  -  v|)  comes  about  due  to  plane  strain  consideratons  (see  section  titled 
“Pure  Bending  of  Moments  in  Plane  Strain”  in  Appendix  C). 


Moment  Computations  for  Bending  Members 

42.  While  use  of  solid  elements  for  bending  members  works  well  to  represent  the 
stiffness  provided  by  bending,  the  problem  remains  as  to  how  to  compute  moments. 

The  solid  element  representation  naturally  provides  statically  equivalent  stress  values  at 
the  center  of  the  element;  these  values  cannot  be  related  to  a  bending  moment.  One 


Figure  5.  Strain  gage  method  of  computing  bending  moments 
for  four-node  solid  element  (Leavell  et  al.  1989) 
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method  used  for  computing  moments  is  based  on  the  premise  that  moments  could  be 
computed  from  beam  theory  using  the  “outer  fiber  strains”  computed  from  displace¬ 
ments  of  the  end  nodes.  This  process  is  illustrated  in  Figure  5,  which  shows  the  solid 
elements  in  a  bending  pattern.  The  outer  fiber  strains  are  shown  to  be  related  to  a 
radius  of  curvature  that  a  true  beam  element  would  conform  to.  As'an  expedient,  the 
outer  fiber  strains  are  computed  by  placing  bar  elements  on  n.o  edges  of  the  beam  ele¬ 
ments.  These  “strain  gage”  elements  are  created  by  using  the  standard  bar  element 
provided  by  SOILSTRUCT  (for  modeling  anchors  and  struts,  etc.).  The  bar  was  given 
a  low  stiffness  so  that  there  was  virtually  no  interaction  between  the  bar  element  and 
surrounding  elements.  The  strains  measured  in  the  two  bars  are  therefore  the  outer 
fiber  strains  £f  and  Ej .  These  strains  may  be  related  to  the  bending  strain  and 
axial  strain  £„  as  follows: 


i<£r 


ei> 


(29) 


£<Cr-  81> 


(30) 


For  the  case  of  pure  bending  (no  axial  load)  £f  =  -  £j  and  £„  *  0  .  For  purely  .xial 
loads  £f»  £j  and  ^  =  0.*  Once  the  strains  have  been  computed  the  moment  j  or 
unit  width  of  sheet-pile  wall  is  obtained  from  the  following  (see  section  titled  "Pure 
Bending  of  Members  in  Plane  Strain”  in  Appendix  C): 


M 


h  (1 


M_e 

-vj)'” 


(31) 


The  factor  of  2  in  the  above  equation  results  from  the  depth  to  neutral  axis  of  nne  half 
the  width,  hfi ,  of  the  corresponding  sheet-pile  element. 

Accuracy  of  Computed  Moments 

43.  Leavell  et  al.  (1989)  investigated  the  ability  of  the  strain-gage  method  to  accurately 
predict  moments  by  comparing  the  momeiits  computed  in  a  finite  element  analysis  of  a 
fixed  end  beam  with  the  moments  computed  using  classical  beam  theory  for  the  same 
cantilever  beam.  They  found  that  when  the  Poisson’s  ratio  is  set  equal  to  0  in  the  finite 


Note  that  a  stiffness  could  be  given  to  the  bar  to  customize  the  beam  element  for 
unsym metrically  reinforced  concrete  walls,  etc.,  or  to  model  tensile  cracking  of  walls  by 
using  a  compression-only  bar.  Also,  pure  shear  deformation  of  the  pile  causes  no  strain  in 
the  bars,  a  fact  that  could  be  of  some  importance  since  the  moment  of  inertia  (I)  scales  as 
the  cube  of  the  pile  thickness  whereas  the  shear  stiffness  is  proportional  to  thickness. 
Thus,  the  bars  could  be  used  to  add  stiffness  to  bending  without  changing  shear  behavior. 


/ 

./ 

/■' 
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element  analysis,  the  computed  results  differed  from  those  results  computed  using  clas¬ 
sical  beam  theory  by  0.01  percent. 

44.  The  displacement  along  the  beam  is  approximated  by  the  solid  element  as  a 
series  of  straight  lines.  (If,  instead,  the  beam  is  represented  by  a  true  bending  element, 
the  displacement  would  be  represented  by  a  smooth  curve.)  As  a  result,  the  bending  mo' 
ment  computed  for  the  element  represents  an  average  value  that  is  presumably  indica¬ 
tive  of  the  value  at  the  center  of  the  element.  Mosher  and  Knowles  (1990)  showed  that 
the  resolution  can  be  improved  by  using  more  elements  to  represent  the  pile  and  by 
using  elements  with  smaller  aspect  ratios  to  model  the  pile.  However,  the  important 
feature  of  the  solid  elements  is  that,  in  contrast  to  standard  beam  elements,  they  deform 
in  a  manner  that  is  compatible  with  adjacent  soil  elements,  a  consideration  of  far 
greater  importance  than  the  small  error  inherent  with  the  linear  approximation. 

Elements  for  Modeling  Sheet  Pile  Cells 

45.  SODLSTRUCT  can  be  used  for  analysis  of  rows  of  sheet  pile  cells  by  considering 
an  “equivalent”  planar  slice  through  the  row  (Clough  and  Kuppusamy  1985).  Deter¬ 
mination  of  properties  for  the  sheet-pile  elements  is  complicated  by  two  aspects  of  the 
problem,  (a)  the  arcuate-shaped  cells  are  modeled  by  a  planar  slice  that  represents  the 
average  behavior  of  the  main  and  arc  cells  and  (b)  the  sheet-pile  properties  represent  an 
average  behavior  of  the  web  and  interlocks. 

46.  To  address  the  first  problem  created  by  the  planar  representation,  the  action  of 
the  sheet  piles  was  depicted  by  Clough  and  Kuppusamy  (1985)  as  consisting  of  three 
independent  parts;  the  outermost  (riverward  and  landward)  sheet  piles  represented  as 
sheet-pile  walls  using  solid  elements,  the  hoop  stiffness  of  the  cell  represented  by 
horizontal  springs  connecting  the  outside  walls,  and  sheet  piles  separating  the  main  and 
arc  cells  represented  as  shear  walls. 

47.  The  justification  for  the  planar  representation  depends  on  the  independence  of 
the  three  cell  components  described  above.  The  independence  of  the  three  components, 
in  turn,  depends  on  the  small  flexural  rigidity  of  the  sheet  piles  that  make  up  the  cells. 
Because  very  little  resistance  is  derived  through  bending,  the  axial  stiffness  of  the  sheet 
piles  is  the  principal  structural  component  of  the  outside  walls.  Therefore,  the  vertical 
resistance  derived  from  the  arcuate-shaped  wall  can  be  modeled  by  a  planar  wall 
reasonably  well.  To  illustrate  how  shear  resistance  is  derived  from  axial  sheet  pile  stiff¬ 
ness,  the  equilibrium  of  the  sheet  pile-soil  system  is  shown  in  Figure  6.  The  ability  of 
the  pile  to  sustain  axial  loads  permits  a  jump  in  shear  stress  across  the  pile  equal  to  the 
gradient  in  axial  load  along  the  pile.  Without  the  pile,  a  continuous  horizontal  shear- 
rupture  plane  could  form,  causing  a  sliding  failure.  With  the  pile  in  place,  the  shearing 
stress  is  distributed  along  the  pile,  and  for  the  cell  to  fail,  vertical  shear  failure  must 
occur  along  the  length  of  the  wall.  The  correspondence  between  the  jump  in  shear 
stress  across  the  sheet-pile  wall  and  the  gradient  of  axial  pile  force  from  a  typical 
analysis  can  be  seen  clearly  in  Figure  6. 
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Figure  6.  Example  of  stress  transfer  from  soil  to  sheet  pile 
(Peters,  Holmes,  and  Leavell  1986) 


48.  The  hoop  stresses  carried  by  the  springs  respond  primarily  to  changes  in  lateral 
loads  and  are  thus  most  important  during  cell  filling,  excavation,  and  flood  loading. 

The  lateral  stress  developed  in  the  cell  fill  and  the  sheet-pile  interlock  forces  depend 
directly  on  the  hoop  stiffness.  The  springs  closely  approximate  the  hoop  stresses, 
provided  the  shear  stress  carried  by  the  sheet  piles  is  small.  It  has  Been  observed  from 
analysis  that  the  spring  stiffness  does  not  directly  influence  the  shear  resistance  of  the 
cell.  The  shear  resistance  of  the  cell  is  derived  from  the  outside  walls  and  the  shear 
resistance  of  the  soil.  The  springs  contribute  to  shear  resistance  by  laterally  confining 
the  cell  fill. 

Outside  wall  solid  elements 

49.  The  properties  of  the  solid  elements  making  up  the  outside  walls  are  selected  to 
match  the  axial  stiffness  for  a  unit  length  of  cofferdam.  For  example,  Peters  et  al. 
(1986)  used  a  solid  sheet-pile  element  1  ft  depth  with  a  modulus  of  4.6  x  10s  psf,  a 
combination  of  thickness  and  modulus  that  gave  an  equivalent  stiffness  of  the  sheet 
piling  while  not  creating  an  excessive  length-to-depth  ratio  for  the  element.  The 
modulus  was  found  from  scaling  the  element  modulus  for  the  planar  model  to  the  actual 
main-arc  cell  combination  based  on  equivalence  of  steel  area.  The  combined  main  cell 
and  arc  consisted  of  approximately  170  PS -32  sheet  piles  with  each  pile  having  an  area 
of  1 1.8  in2.  The  total  pile  area  is  2,006  in2  or  about  14  ft2.  The  length  represented  by 
an  arc  and  main  c  d  is  approximately  65  ft,  giving  a  unit  steel  area  of  0.22  ft2/ft.  The 
corresponding  unit  steel  area  for  the  two  1-ft-thick  elements  is  2  ft2/ft.  To  compensate 
for  the  greater  area  in  the  planar  model,  the  modulus  must  be  multiplied  by  a  factor  of 
0.11.  Assuming  a  modulus  for  steel  of  29  x  106  psi,  the  element  modulus  should  be 
3.19  x  106  psi,  or  approximately  4.6  x  108  psf. 

Hoop  stiffness  spring  elements 

50.  The  spring  stiffness  is  determined  by  assigning  values  that  would  make  the 
response  of  the  planar  cell  to  an  internal  radial  pressure  equivalent  to  the  response  of  a 
circular  cell.  If  a  circular  cell  is  subjected  to  a  uniform  internal  pressure,  its  radial  dis¬ 
placement  u  is  given  by: 


u 


P  x  R2 
E  x  t 


(32) 


where 

Pa  radial  pressure  inside  cell 
R  »  cell  radius 
t  a  cell-wall  thickness 
which  gives  an  equivalent  spring  stiffness  of: 


p 

u 


Et 


(33) 


K  = 


where  K  represents  a  uniformly  distributed  radial  spring  acting  to  restrain  radial  move¬ 
ment.  The  planar  structure  is  to  be  represented  by  two  cell  walls  connected  by  discrete 
springs.  The  total  effect  of  the  springs  in  resisting  lateral  movement  over  a  vertical 
tance  H  is: 

Ks-  4-  (34) 

*  6 


where  6  is  the  lateral  movement  of  the  wall  and  p  is  the  resultant  force  PH. 

The  relationship  between  the  radial  displacement  of  the  circular  cell  and  the  lateral 
movement  of  the  two  walls  is  5  =  2u  .  Thus,  to  obtain  the  equivalent  displacement  for 
a  given  pressure: 


K 


s 


(35) 


For  the  example  considered  by  Peters  et  al.  (1986),  the  stiffness  K  can  be  determined 
from  the  sheet-pile  properties;  E  =  29  x  106  psi ,  v  =  0.3  ,  and  t  =  0.5  in. ,  and  cell 
diameter  2R  =  52  ft ,  to  get 

K  =  3.03  x  105  pcf 
which  gives  a  planar  spring  stiffness  of: 

Ks  =  1.51  x  105  •  H  pcf  (3 


where  H  is  the  vertical  spring  spacing. 

51.  The  spring  constant  derived  above  was  based  on  the  assumption  that  the  effec¬ 
tive  modulus  of  the  sheet  piling  was  equal  to  that  of  a  steel  membrane  having  a  thick¬ 
ness  t .  In  fact,  displacement  of  the  interlocks  reduces  the  modulus  in  the  tangential 
(horizontal)  direction.  To  account  for  interlock  behavior,  the  modulus  is  multiplied 
by  an  “E  ratio”  which  is  defined  as  the  ratio  of  the  effective  modulus  to  the  modulus 
of  solid  steel.  An  E  ratio  of  0.03  has  been  found  from  experience  (Clough  and 
Kuppusamy  1985)  to  yield  a  reasonably  good  correspondence  between  model  and 
field  behavior.  Thus  the  spring  modulus  used  in  the  above  example  should  be: 
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(37) 


K.  =  4.5  x  103  •  H  pcf 
s 


In  addition  to  determining  a  conect  spring  constant,  it  is  also  necessary  to  relate  the 
spring  force  to  the  interlock  stress.  The  strain  in  the  circular  cell  is  given  by: 

_  _  _u  _  P  x  R  f.  _  ij]  (38) 

"  R  "  E  x  t  [  2 


Ignoring  the  strain  in  the  axial  direction,  the  horizontal  stress  is  given  by  c  =  Ee .  The 
horizontal  force  per  unit  length  of  interlock  is  T  =  crt .  The  unit  interlock  force  is 
therefore  given  by: 


=  0.85  P  x  R 


(40) 


Following  a  similar  line  of  reasoning,  the  spring  force  T$  can  be  related  to  the  displace¬ 
ment  by  Ts  =  2Ksu  and  to  the  interlock  force  by  T  =  T,  /H  .  Therefore  by  substitut¬ 
ing  K  =  Ks /H  ,  the  equivalence  between  T  and  Ts  can  be  developed  as  follows: 
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Figure  7.  Common-wall  model  ( Peters ,  Holmes,  andLeavell  1986) 
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Thus,  in  this  example  the  unit  interlock  force  is  related  to  the  spring  force  by: 

T  R 

T  =  0.85  -2—  lb/ft 

H  * 

Alternative  description  of  common-wall  elements 

52.  The  common  wall  sheet  piles  separating  the  main  and  arc  cells  provide  shear 
resistance  that  is  not  conveniently  accounted  for  when  modeled  using  spring  elements. 
To  better  approximate  shear  stiffness  of  the  common  wall,  the  soil  elements  can  be  over- 
lain  with  “patch  elements”  which  act  independently  of  the  soil  except  where  attached  at 
the  outside  walls  (Peters,  Holmes,  and  Leavell  1986).  Figure  7  gives  a  graphical 
description  of  how  this  technique  is  implemented.  The  bar  elements  account  for  the 
higher  proportion  of  steel  area  in  the  outside  wall,  providing  additional  axial  stiffness 
without  contributing  to  flexural  or  shear  stiffness.  The  interface  elements  model  the 
friction  between  the  soil  and  outside  sheet-pile  wall,  as  in  the  spring  model. 

53.  The  patch  elements  are  characterized  as  a  transversely  isotropic  material  acting 
under  plane  stress  conditions.  Five  material  properties  are  required  to  define  their  be¬ 
havior.  These  properties  are:  Ej  and  Vj  ,  which  are  associated  with  the  common- 
wall’s  axial  behavior  and  E2  ,  v2 .  and  the  shear  modulus  G2 ,  which  define  stiffness 
in  the  tangential  direction.  The  subscripts  follow  the  convention  used  by  Zienkiewicz 
(1977).  These  properties  allow  not  only  the  shear  stiffness  to  be  modeled  but  also  the 
difference  between  the  vertical  and  horizontal  stiffness  of  the  cell.  Note  that  Vj  does 
not  appear  in  the  plane  stress  stiffness  formulation  and  v2  can  be  set  to  equal  zero  to 
simplify  property  determination.  In  fact,  the  interlocks  may  tend  to  inhibit  the  Poisson’s 
effect,  although  no  data  exist  to  support  this  assumption.  Provisions  have  been  made  in  the 
finite  element  code  for  and  values  other  than  zero  could  be  included. 

Development  of  anisotropic  properties 

34.  The  rationale  for  estimating  the  stiffness  properties  follows  the  same  line  of 
reasoning  outlined  previously:  the  various  components  of  stiffness  can  be  uncoupled 
and  treated  independently.  The  reasoning  for  each  parameter  is  described  as  follows: 

55.  E2.  The  equation  for  determining  the  hoop  stiffness  E2  is  derived  by  relating 
the  response  of  a  vertical  planar  slice  through  the  center  of  the  cell  to  that  of  a  circular 
cell.  Subjecting  a  circular  cell  to  a  uniform  internal  pressure,  its  radial  displacement  u 
is  given  by: 


where 
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E,  =  Young’s  modulus  of  steel 
v  =  Poisson’s  ratio,  assumed  to  equal  0.0 

The  lateral  displacement  8  caused  by  the  same  pressure  in  the  planar  slice  is  given  by: 

8  =  Jj 1  _L  (*?> 


where 


P  =  lateral  pressure 
d  =  unit  thickness  of  slice 

Ej  =  hoop  stiffness  or,  equivalently,  the  effective  modulus  of  the  slice 
t  =  slice  thickness  (=  1) 

L  s  slice  length 

Substituting  5  =  2u  ,  L  =  2R  ,  and  P  =  P  into  Equation  47  and  equating  Equations 
46  and  47,  it  is  found  that 


Esx  t 

R 


(48) 


Solving  Equation  48  with  the  following  variable  values  (Es  =  4.18  x  109  psf , 
t  =  0.0667  ft ,  and  R  =  25  ft)  results  in  a  hoop  stiffness  of  E2  =  1.12  x  107  psf  . 
Because  of  displacement  in  the  sheet-pile  interlocks,  E2  is  multiplied  by  an  “E  ratio” 
of  0.03,  giving 

E2  =  3.34  x  105  psf 

56.  The  procedure  outlined  above  is  clearly  equivalent  to  that  used  to  determine  the 
horizontal  spring  stiffness  in  the  conventional  planar  model.  If  Ej  and  G2  are  set  to 
zero,  the  patch  elements  behave  as  a  continuous  spring. 

57.  E|  and  G2  .  The  vertical  stiffness  Ej  is  obtained  by  determining  the  percent 
steel  in  a  typical  planar  section  of  the  sheet  pile  wall  and  adjusting  the  modulus  of  steel 
to  obtain  an  equivalent  stiffness.  Using  the  total  length  of  sheet  pile  material  as  240  ft, 
an  average  thickness  of  0.0667  ft,  and  a  total  area  of  3,176  ft3,  the  percent  steel  in  a  sec¬ 
tion  of  wall  is  approximately  0.51  percent.  Multiplying  the  modulus  of  steel,  4.18  x 
)09  psf,  times  0.51  percent  gives 

Ej  =  2.13  x  107  psf 
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The  shear  modulus  G2  is  assumed  to  be  0.4  of  Ej ,  making 
G2  =  8.52  x  106psf 

The  0.4  multiplier  accounts  for  the  Poisson's  ratio  of  the  steel  within  the  pile  web.  The 
present  model  does  not  account  for  slip  along  the  interlocks  which  would  tend  to  reduce 
the  shear  stiffness. 

Axial  bar  stiffness 

58.  The  axial  bar  placed  along  the  outer  walls  accounts  for  the  fact  that  the  area  of 
steel  at  those  locations  is  the  same  in  both  the  prototype  and  the  planar  slice.  The  area 
of  steel  computations  for  the  bar  stiffness  is  identical  to  that  described  for  the  solid  wall 
elements  in  previous  sections  of  the  report. 


PART  IV:  USE  OF  SOILSTRUCT  PROGRAM 
Sign  Convention  and  Coordinate  System 

59.  All  input  data  and  results  are  specified  using  a  right-hand  coordinate  system;  the 
x-axis  being  horizontal  and  positive  to  the  right  and  the  y-axis  being  vertical  and  posi¬ 
tive  upwards.  The  sign  convention  for  stresses  acting  at  the  center  of  a  two-dimen¬ 
sional  element  is  shown  in  Figure  8a.  Compressive  stresses  are  taken  to  be  positive. 

60.  Stresses  for  interface  elements  are  defined  with  respect  to  their  local  axes  along 
the  length  of  the  interface  x'  as  defined  by  the  I  and  J  nodes,  and  normal  to  the  ele¬ 
ment  /  as  shown  in  Figure  8b.  Positive  normal  stresses  are  compressive.  Positive 
shear  stresses  act  in  the  positive  x'  direction  along  the  length  of  the  interface  as  shown 
in  Figure  8b. 

61.  Positive  forces  are  taken  to  be  compressive  in  all  one-dimensional  bar  elements, 
with  the  exception  of  bar  reinforcement  elements.  For  these  elements,  positive  forces 
are  taken  to  be  tensile. 

Units 

62.  Any  consistent  set  of  units  can  be  used  with  SOILSTRUCT,  with  one  cautionary 
note.  The  normal  stiffness  of  interface  elements  is  arbitrarily  set  to  a  value  of  1  x  108, 
independent  of  units,  as  discussed  in  paragraph  33 . 

Capacity 

63.  The  capacity  of  the  program  is  determined  by  the  size  of  variables  in  the  com¬ 
mon  block  and  dimension  statements.  The  global  stiffness  matrix  [K]  is  stored  within 
the  program  in  vector  SN  .  The  dimension  of  SN  is  equal  to  1,600,000,  and  the  maxi¬ 
mum  number  of  either  nodes  or  two-dimensional  elements  is  4,000.  The  dimensions  for 
all  of  the  arrays  are  given  throughout  the  input  guide.  Double  precision  is  specified  for 
all  arrays  in  the  CRAY  Y-MP  version  of  SOILSTRUCT. 

64.  The  global  stiffness  matrix  is  assembled  in  vector  form  using  a  modified  skyline 
procedure.  Therefore,  there  is  no  direct  correlation  between  the  size  of  SN  and  the  max¬ 
imum  number  of  nodes  or  elements.  If  the  required  size  of  SN  exceeds  1,600,000 
during  execution  of  SOILSTRUCT,  the  program  will  terminate  execution  and  print  the 
size  required  for  the  SN  vector.  If  this  occurs,  two  statements  must  be  modified 
within  the  computer  code;  (a)  the  dimension  of  the  SN  vector,  found  in  the  main  pro¬ 
gram,  must  be  increased  to  the  required  value,  and  (b)  a  call  statement,  comparing  the 
required  size  of  the  SN  vector  for  the  problem  being  analyzed  to  the  actual  dimension 
of  SN ,  must  be  revised  in  subroutine  DETNA. 
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°y 


a.  Coordinate  system  and  sign  convention  for  positive  stresses 
for  two-dimensional  elements 


y.x' 


L  I  STRESS-7 

b.  Definition  of  interface  local  axes  and  sign 
convention  for  positive  shear  stress 


Figure  8.  Sign  conventions  and  coordinate  axes 
for  stress  and  displacements  (Clough  1984) 
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APPENDIX  A:  USER’S  GUIDE  FOR  PROGRAM  SOILSTRUCT 


Input  Data  Sequence 


The  format  of  the  input  data  is  free-field  format  unless  stated  otherwise. 
1.  IDENTIFICATION  CARD  —  FORMAT  (20A4) 


Column  Variable 

1-80  HED 


Explanation 


Analysis  identification. 


2a.  DATA  CONTROL  CARD 


All  nodes  and  two-dimensional  and  interface  elements  to  be  used  in  the  analysis  must  in¬ 
itially  be  included  in  the  mesh;  additions  or  deletions  of  these  elements  are  not  allowed, 
but  the  material  parameters  can  be  changed  to  make  them  inactive.  One-dimensional  bar 
elements  may  initially  be  included  in  the  mesh  or  added  in  subsequent  construction  steps. 


Variable 

NUMNP 

NUMEL 


NUMJT 

NUMBAR 


NC 

NMOD 


Number  of  nodal  points  (4,000  maximum). 

Number  of  elements  initially  in  the  mesh,  exclud¬ 
ing  bar  and  beam  elements,  but  including  interface 
elements.  NUMEL,  then,  includes  NUMJT  (4,000 
maximum).  Interface  elements  should  be  num¬ 
bered  first. 

Number  of  interface  elements  (400  maximum). 

Number  of  bar  elements  initially  in  the  mesh,  in¬ 
cluding  those  bar  elements  added  in  subsequent 
construction  stages  (400  maximum).  NUMBAR 
also  includes  NSTRAN,  the  number  of  bar  elements 
used  as  strain  gages  (see  card  number  2c). 

Number  of  loading  and  construction  steps  (40 
maximum). 

Modulus  specification  code- 

■  0  if  modulus  calculation  codes  input  with 

loading  information  card. 

■  1  if  modulus  calculation  codes  input  with 

modulus  calculation  card. 
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INIT  Initial  stress  input  code  - 

=  0  if  external  input  from  cards  or  tape,  included 
in  input. 

=  1  if  internally  generated  from  gravity  turn-on 
analysis. 

=  2  if  initial  stresses  and  displacements  are  to  be 
set  equal  to  zero  within  SOILSTRUCT. 


KI 


If  INIT  *  0  or  *  1  or  *  2  ,  initial  stresses  are  generated  assuming  a 
horizontal  ground  surface,  horizontal  water  table,  and  ax  =  KQ  oy  . 


Interface  element  activation  code  - 
=  0  not  activated  during  initial  stress  computation. 
=  1  activated  during  initial  stress  computation. 

KI  is  used  when  INIT  =  0  .  If  INIT  =  0  ,  then  KI 


can  be  set  equal  to  zero. 


IHORIZ 


ITRD 


Ground  surface  inclination  code  - 
=  0  horizontal  ground  surface.  Vertical  stresses 
are  computed  using  a  gravity  turn-on  method 
of  analysis.  Horizontal  stresses  are  computed 
assuming  K0  =  v/(l-  v)  unless  KQ  is  specified. 
=  1  sloping  ground  surface.  Vertical  and  horizontal 
stresses  are  calculated  from  a  gravity  turn-on 
analysis  assuming  linear  elastic  response  of  soil 
-  i.|e.,  K0  =  v/(l-  v) . 

Analysis  printout  code  - 
=  0  if  initial  stresses  and  results  of  the  final 
iteration  are  to  be  printed. 

=  1  if  initial  stresses  and  results  of  all  iterations 


are  to  be  printed. 

=  2  if  initial  stresses  are  not  printed  but  results  of 
all  iterations  are  to  be  printed. 

=  3  if  initial  stresses  and  results  of  the  final 


iteration  are  to  be  printed.  No  substep  results 
are  to  be  printed. 

=  -2  if  initial  stresses  are  not  printed  but  results  of 
the  final  iteration  are  to  be  printed. 


ILIST  Element  and  nodal  point  card  data  printout  code 

=  0  if  not  printed. 

=  1  if  printed. 
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Variable 


Explanation 


DPUNCH  Restart  file  code - 

=  0  if  no  restart  file  is  to  be  generated. 

*  1  if  a  restart  file  is  to  be  generated.  Only  the 
results  of  the  final  iteration  are  used. 

IT  APE  Disk  storage  code  - 

=  0  if  no  disk  storage. 

**  1  if  storage  of  displacements  from  final  iteration. 
=  2  if  storage  of  all  data  except  displacements 
from  final  iteration. 

GAMW  Unit  weight  of  water. 

PATM  Atmospheric  pressure. 

The  unit  weight  of  water  and  the  atmospheric  pressure  are  included  as  basic  parameters. 
Either  English  or  SI  units  can  be  used.  All  data  must  be  compatible  with  input  coor¬ 
dinate,  pressure,  and  material  property  parameters. 

2b.  MATERIAL  ALLOCATION  CARD 

All  two-dimensional  material  types  are  assigned  material  numbers  first,  followed  by 
the  interface  material  types.  Bar  elements  are  not  assigned  a  material  type  number,  but 
are  identified  solely  by  their  element  number.  If  a  number  for  NATYP,  or  NA2TYP, 
for  example,  is  not  required,  assign  a  value  of  zero. 


Variable 

Explanation 

NUMMAT 

Total  number  of  material  types,  including  both 
two-dimensional  soil  or  construction  material 
types  and  interface  material  types. 

NUMSOL 

Total  number  of  material  types  excluding  the  inter¬ 
face  material  types.  Thus,  (NUMMAT-NUMSOL) 
must  equal  the  number  of  interface  material  types. 

NATYP 

Material  type  number  assigned  to  elements  initially 
having  the  properties  of  air.  Usually,  elements  that 
will  be  built  are  initially  identified  as  air  elements. 

NCTYP 

Structural  material  type,  such  as  concrete  or  sheet¬ 
piling. 

NB1TYP 

Backfill  material  type  1.  (Refer  to  section  13b  on 
fill  or  concrete  placement.) 

NB2TYP 

Backfill  material  type  2.  (Refer  to  section  13b  on 
fill  or  concrete  placement.) 

Variable 


Explanation 


NSTRAN 

NOMOM 

EBAR* 

EBEND* 

NA2TYP 

R_FILL 


Number  of  bar  elements  used  as  strain  gages. 
These  bar  elements  are  numbered  first. 

Number  of  pairs  of  bar  element  strain  gages  used 
in  moment  computations  (see  paragraph  42). 

Young’s  modulus  for  steel  reinforcement.  Used  in 
conjunction  with  bar  element  strain  gages. 

Young’s  modulus  for  the  bending  member.  Used 
in  conjunction  with  bar  element  strain  gages  and 
moment  computations. 

Material  type  number  assigned  to  interface  ele¬ 
ments  having  the  properties  of  air. 

Total  number  of  soil  element  material  types  in 
which  layers  of  reinforcement  are  present.  When 
conducting  a  restart  analysis  in  which  reinforce¬ 
ment  is  present,  specify  a  negative  value  for 
R_FILL. 

If  any  of  the  above  are  not  needed,  enter  zero. 


°  Set  equal  to  Ee  for  sheet  piles  (see  paragraph  41). 
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2c.  CALCULATION  PROCEDURES 


These  five  integers  activate  special  features  incorporated  within  this  version  of 
SOILSTRUCT. 

Variable  _ Explanation _ 


KEYEI 


Procedure  used  to  define  the  parameter  MQ  and 
the  shear  modulus  MD  of  the  stress-strain  matrix 
[D]  where  {o}  =  [D]{e} . 


(MB  +  Mjj)  (MB  -  Md)  0 
(Mg  -  Mjj)  (Mg  +  Md)  0 
0  0  Md 


=  -1,  then 

Mg  =  Mg(Ej,v) 
MD  =  MD(Et,v) 


=  0, then 

Mg  =  Mg(Et,B) 
Md  =  Md(E,,B) 


=  l.then 

Mg  *  Mb(Ej,B) 
Md  =  MD(Et,B) 


IPETER  Procedure  used  to  calculate  the  stress  level,  SL 

(see  Figure  2). 

=  0  when  SL  is  calculated  using 
(o3>failure  =  current  ' 

=  1  when  SL  is  calculated  using  mean  pressure  at 
failure  equal  to  current  mean  pressure. 

NPETER  Procedure  used  to  update  the  total  interface  nodal 

point  displacements. 

■  0  if  I,  J,  K,  and  L  interface  node  displacements 
are  updated. 

=  1  if  interface  node  displacements  I  and  L  are 
set  equal,  and  node  displacements  J  and  K 
are  set  equal.  This  option  is  used  for  new  inter¬ 
face  elements  placed  adjacent  to  pre-existing 
structural  elements. 


AS 


KEYH20  Procedure  used  to  define  the  lateral  stress  within 

newly  placed  fill  elements  (see  card  group  13b). 

=  Owhen  cx=  K0<jy 

=  1  when  cr  =  ov 

A  y  • 

XEYSEP  Applied  seepage  forces  (see  card  group  13c). 

=  0  if  lateral  and  vertical  seepage  forces  are  to  be 
applied. 

=  1  if  only  vertical  seepage  forces  are  to  be  applied. 
3.  LOADING  INFORMATION  CARD 


One  card  is  supplied  for  each  load  step.  One  to  three  load/  construction  modes  can 
be  included  in  each  load  step.  The  load  or  construction  mode  codes  include: 

KCS(NCJ)  _ DESCRIPTION _ 

1  Excavation  (eauivalent  nodal  loads  can  be  applied 
in  equal  increments). 

2  Fill  placement  (subroutine  SUBSTP  cannot  be 
used  in  conjunction  with  the  fill  placement  proce¬ 
dures  of  subroutine  BUILD). 

3  Seepage  loading  (equivalent  nodal  loads  can  be  ap¬ 
plied  in  equal  increments). 

4  Deletion  of  bar  element  (force  in  the  element  can 
be  applied  in  equal  increments). 

5  Installation  of  bar  element  (prestress  force  can  be 
applied  in  equal  increments). 

6  Boundary  pressure  loading  (equivalent  nodal  loads 
can  be  applied  in  equal  increments). 


7  Temperature  loading  (the  total  temperature  change 
can  be  applied  in  equal  increments).  1 

8  Support  displacement  (the  total  displacement  can 
be  applied  in  in  equal  increments). 

9  Concentrated  nodal  loads  (can  be  applied  in  equal 

increments).  \ 

10  Element  material  type  change.  1 

As  indicated  in  the  listing,  prescribed  loads,  displacements,  or  temperature  changes  can 
be  analyzed  in  equal  increments,  or  substeps,  for  each  load  case.  Subroutine  SUBSTP 
generates  the  equivalent  load  increments,  then  main  analyzes  all  increments  prior  to 
analyzing  the  next  load  step.  With  one  exception,  all  loading/construction  modes  that 
can  be  applied  in  increments,  or  substeps,  can  also  be  applied  in  any  combination  in  any 
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load  step.  The  number  of  substeps,  however,  will  be  the  same  for  all  loading  or  con¬ 
struction  modes  included  in  the  load  step.  The  exception  is  temperature  loading;  if  a 
temperature  change  is  specified,  and  a  given  number  of  substeps  are  specified,  then 
only  the  corresponding  temperature  loads  can  be  specified  in  the  loading  step  -  i.e., 
KCS(N,2)  and  KCS(N,3)  must  be  set  equal  to  zero.  If  the  number  of  substeps, 
NSBSP ,  is  equal  to  zero,  then  temperature  loading  can  be  included  with  other  loading/ 
construction  modes  in  a  load  step. 

Since  the  same  input  format  is  used  in  modes  8  and  9,  the  following  rules  apply;  if  only 
concentrated  nodal  loads  are  specified,  use  mode  9,  if  only  support  displacements  are 
specified,  use  mode  8,  and  if  both  loads  and  displacements  are  specified,  use  mode  8. 


Variable 


KCS(NC,1)  First  loading/construction  mode  code. 


KCS(NC,2)  Second  Ioading/construction  mode  code. 


KCS(NC,3)  Third  loading/construction  mode  code. 


KSC(NC,1) ,  KCS(NC,2) ,  and  KCS(NC,3)  can 
be  input  in  any  numerical  order,  but  the  modes  are 
processed  in  ascending  numerical  order:  If  the 
second  and/or  third  loadings  are  not  required,  then 
KCS(NC,2)  KCS(NC,3)  should  be  set  equal  to 
zero. 


NUMIT(NC)  Number  of  iterations  for  the  load  step. 

NUMIT(NC)  applies  to  each  substep  if  substeps 
are  specified.  NUMIT(NC)  =  0  is  the  same  as 
NUMIT(NC)  =  1 . 

NUMSS(NC)  Number  of  substeps,  the  minimum  value  which 

may  be  assigned  to  NUMSS(NC)  is  one. 

MOD(l,NC)  Modulus  calculation  code- 

=  1  if  a  loading  modulus  is  to  be  calculated. 

=  2  if  an  unload-reload  modulus  is  to  be 
calculated. 

=  0  if  the  computer  is  to  decide  the  type  of 
modulus  to  be  calculated.  In  this  case,  if  the 
most  recently  calculated  maximum  shear  stress 
for  an  element  is  less  than  all  previous  values 
of  maximum  shear  stress,  an  unload-reload 
modulus  is  assigned.  Otherwise  a  primary 
loading  modulus  is  assigned. 


Additional  input  is  required  here  only  if  NMOD  =  0 .  All  material 
types,  other  than  interface  or  bar  elements,  are  given  one  of  the 
above  codes.  If  NMOD  =  0  and  NC  =  0 ,  as  might  be  the  case 
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for  an  analysis  of  initial  stresses,  M0D(I,1)  is  set  equal  to  zero, 
or  the  computer  decides. 

IPRT(NC)  =  0  Do  not  print  the  force  vector. 

=  1  Print  the  force  vector.  . 

IPLT(NC)  -  0  Do  not  create  a  plot  file. 

=  1  Create  a  plot  file,  including  moments. 

=  2  Create  a  plot  file  of  moments  only. 

=  3  Create  a  NISA  plot  file  of  geometry,  displace¬ 
ments,  and  stresses. 

HEDCS(NC)  Description  of  the  load  step. 

4.  MODULUS  SPECIFICATION  CARD 


This  card  is  required  only  if  NMOD  =  1  and  NC  =  0 .  A  card  is  required  for  each 
loading  step,  1  to  NC  .  In  this  option,  values  of  the  modulus  specification  code  are 
specified  for  each  material  type  (and  thus  each  element,  excluding  bar  and  interface 
elements),  regardless  of  the  change  in  maximum  shear  strain  that  may  have  occurred. 

Variable  _  Explanation  _ _ 

MOD(I,NC)  Modulus  calculation  code  for  each  material  type 

(1  to  NUMSOL)  for  the  first  load  step.  Separate 
cards  are  required  for  each  load  step. 

5.  MATERIAL  PROPERTY  CARDS 

These  pairs  of  cards  are  used  only  to  define  the  material  properties  for  two- 
dimensional  elements,  excluding  bar  and  interface  material  types.  The  first  and  second 
lines  are  supplied  in  order  of  material  type  number  N  =  1  to  NUMSOL  .  Information 
or  properties  not  required  for  a  material  type  can  be  set  equal  to  zero. 

Variable  _ Explanation _ 

GUI(N)  Poisson’s  ratio  before  failure,  or  the  bulk  modulus 

number  Kb . 

GUF(N)  Poisson’s  ratio  at  failure  (no  greater  than  0.49),  or 

the  bulk  modulus  exponent,  m . 

GAM(N)  Total  or  buoyant  unit  weight  (always  specified, 

regardless  of  drained  or  undrained  material  be¬ 
havior). 

FR(N)  The  failure  ratio  Rj. 

AO(N)  Coefficient  of  lateral  earth  pressure  at  rest  K0  as 

pertaining  to  effective  stresses. 
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V  - 


Variable 


Explanation 


PHI(N) 

Friction  angle  in  degrees.  If  IDRAIN  =  2 ,  PHI 
should  be  set  equal  to  zero. 

XXP(N) 

The  modulus  exponent  n .  For  a  linear  elastic 
material  n  must  be  set  equal  to  0.  For  saturated 
soils  when  PHI(N)  =  0. ,  n  is  normally  set  equal 
to  0.00001. 

IDRAIN(N) 

Material  behavior  code  - 
*  -1  if  patch  element 
a  0  if  undrained 
«  1  if  drained 

a  2  if  value  entered  for  cohesion  is  equal  to  the 
ratio  (Su/P;  ) , 

where  p'  =  1/2  (a'lc  +  a'3  ^  . 

HCOEF(N) 

The  modulus  number  K. 

ULCOEF(N) 

The  unload-reload  modulus  number  Kur . 

COHE(N) 

Undrained  shear  strength  or  cohesion.  When 

I  DRAIN(N)  a  2,  COHE(N)  is  set  equal  to  the 
normalized  strength  ratio  (S^/p'  ) . 

E(N) 

Tangent  modulus  at  failure  for  isotropic,  nonlinear 
materials,  or  Young’s  modulus  for  elastic 
materials. 

ALPHA(N) 

Coefficient  of  linear  thermal  expansion  for  struc¬ 
tural  element.  For  non-structural  material  types, 
set  ALPHA(N)  equal  to  zero. 

EIMN(N) 

Minimum  initial  tangent  modulus  for  isotropic, 
nonlinear  materials.  Set  EIMN(N)  equal  to  zero 
for  elastic  materials. 

TENS(N) 

Minimum  allowable  value  of  the  minor  principal 
stress  for  isotropic,  nonlinear  materials.  If  tensile, 
input  TENS(N)  as  a  negative  value.  For  elastic 
materials,  set  equal  to  zero. 

The  variables  described  above  correspond  to  typical  soil  materials.  However,  accord¬ 
ing  to  the  value  given  to  IDRAIN  these  parameters  may  take  on  different  meanings. ' 
The  table  below  gives  the  interpretation  given  to  each  variable  for  each  value  of 
IDRAIN .  Note  that  *  implies  that  the  variable  is  not  used  but  a  value  of  0.0  should  be 
inserted  as  a  “place  holder”  to  insure  other  values  are  read  properly. 
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IDRAIN  = 


IDRAIN  = 


Variable  _ Explanation 


0  (Total  stress  “undrained”  specification) 


GUI 

GUI  (.49+) 

GUF 

* 

GAM 

GAM  (total) 

FR 

FR 

AO 

K0  (drained) 

PHI 

PHI  (total) 

XXP 

XXP  (>  lOe  -  5) 

IDRAIN 

0 

HCOEF 

HCOEF 

ULCOEF 

ULCOEF 

COHE 

COHE  (total) 

E 

* 

ALPHA 

ALPHA 

EIMN 

EIMN 

TENS 

TENS 

(Effective  stress  specification) 

GUI 

GUI  or  Kb 

GUF 

GUF  or  m 

GAM 

GAM  (bouyant  if  below  water  table) 

FR 

FR 

AO 

K0  (drained) 

PHI 

PHI  (effective) 

XXP 

XXP 

IDRAIN 

1 

HCOEF 

HCOEF 

ULCOEF 

ULCOEF 

COHE 

COHE  (effective) 

E 

* 
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IDRAIN  = 


IDRAIN  = 


't'  '  - 

Variable 

*  *  H  ,  .  , 

*  *  .  /  !  J  '  " 

Explanation 

ALPHA 

ALPHA 

EIMN 

EIMN 

TENS 

TENS 

(Elastic  materials; 

used  for  MTYPE  =  NCTYP  OR  NATYP ) 

GUI 

GUI  or  Kb 

GUF 

* 

GAM 

GAM 

FR 

* 

AO 

K0  (drained) 

PHI 

* 

XXP 

0 

IDRAIN 

0  or  1 

HCOEF 

* 

ULCOEF 

* 

COHE 

* 

E 

E 

ALPHA 

ALPHA 

EIMN 

* 

TENS 

* 

(Undrained  with  properties  based  on  initial  consolidation  stress) 

GUI 

GUI  (=.49+) 

GUF 

* 

GAM 

GAM  (Total  above  water  table  &  buoyant  below) 

FR 

FR 

AO 

K0  (drained) 

PHI 

PHI  (0.) 

XXP 

XXP  (small  but  >  lOe  -  5) 

IDRAIN 

2 

HCOEF 

E/su 

ULCOEF 

Eu/Su 

All 


Variable 


Explanation 


COHE 

VPc 

E 

* 

ALPHA 

* 

EIMN 

EIMN 

TENS 

TENS/SU 

IDRAIN  =  -1  (Anisotropic 

“patch”  elements  used  to  model  sheet  pile  cells) 

GUI 

v2 

GUF 

* 

GAM 

* 

FR 

* 

AO 

* 

PHI 

Et/Ej 

XXP 

g2/e2 

IDRAIN 

-1 

HCOEF 

* 

COHE 

* 

E 

% 

ALPHA 

* 

EIMN 

* 

TENS 

* 

6.  INTERFACE  PROPERTY  CARD 

One  card  is  supplied  for  each  interface  material  type,  N  =  1  to  NUMJT .  If  no  inter¬ 
face  elements  are  used,  no  cards  are  required.  RKN  (N)  is  set  equal  to  1  x  10®  (or  =  1 
if  air  interface  element)  within  the  code. 

Variable 

Explanation 

PHJ(N) 

Interface  friction  angle  in  degrees. 

RKS(N) 

The  minimum  value  for  the  interface  shear  stiff¬ 
ness  k$. 

COJ(N) 

Interface  cohesion  c . 

FRJ(N) 

Failure  ratio  Rfi . 
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Variable 


Explanation 


TENSJ(N) 

IADJMT(N) 

RKJ(N) 

XXPJ(N) 


Tensile  strength  of  interface  ot . 

Adjacent  two-dimensional  element  material  type 
number. 

The  interface  modulus  number  Kj . 

The  interface  modulus  exponent  n4 . 


7.  NODAL  POINT  CARDS  —  FORMAT  (110,  4D10.4) 

One  card  is  supplied  for  each  node.  The  numbering  of  nodal  points  must  be  sequen¬ 
tial  and  the  input  for  some  of  the  nodes  can  be  omitted.  Those  nodes  omitted  are  auto¬ 
matically  generated  by  the  program  at  equal  spacings  between  those  nodes  specified. 
The  first  and  last  nodes  must  always  be  specified.  Note  that  DP(N)  and  PP(N)  are 
automatically  generated  in  equal  increments  for  those  nodes  omitted  (see  Figure  Al). 


Column 

Variable 

Explanation 

1-10 

N 

Nodal  point  number. 

11-20 

X(N) 

X  coordinate,  positive  to  right. 

21-30 

Y(N) 

Y  coordinate,  positive  upward. 

31-40 

PP(N) 

Pore  pressure  in  head  of  water,  zero  or  blank  if  not 
specified.  Pore  pressure  must  be  specified  for 
drained  materials  but  not  for  undrained  materials. 

41-50 

DP(N) 

Change  in  pore  pressure  in  head  of  water  for  soil 
elements;  change  in  temperature  for  linear  elastic 
structural  material,  otherwise  set  equal  to  zero. 

8.  BOUNDARY  CONDITION  CARDS 


Cards  1  through  8  are  supplied  as  required  to  specify  restraints  of  boundary  nodes.  If 
there  are  no  restrained  nodes,  then  this  series  of  cards  is  not  required.  On  each  card,  the 

_ * _ • _ i _ i _ _ a  i _ •  _ » _ l _ 


constrained  nodes  must  be  in  sequential  order, 
(a)  FIRST  CARD 
Variable 


NOY 


NOX 


Explanation 


lumber  of  nodal  points  fixed  against  Y-movement 
?nly. 

lumber  of  nodal  points  fixed  against  X-movement 
only. 


NOXY 


Number  of  nodal  points  fixed  against  both  X-  and 
Y-movement. 
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PRIMED  NUMBERS 
INDICATE  NODE  OR 
ELEMENT  NUMBERS 
GENERATED  BY  THE 
COMPUTER.  USING 
THE  INPUT  (UNPRIMED) 
NODE  AND  ELEMENT 
NUMBERS. 


Figure  Al.  Examples  of  numbering  nodes  and  elements 
if  generated  by  the  program  ( Clough  1984) 


(b)  SECOND  CARD 
Variable 
IC(N) 


(c)  THIRD  CARD 
Variable 
IC(N) 


(d)  FOURTH  CARD 
Variable 
IC(N) 


9.  ELEMENT  CARD 


Nodal  point  number  of  the  first  nodal  point  fixed 
against  Y-movement.  Additional  nodal  points 
fixed  against  Y-movement,  N  =  2  to  NOY ,  are 
specified  in  the  next  columns  and  on  additional 
cards  as  required. 


Nodal  point  number  of  the  first  nodal  point  fixed 
against  X-movement.  Additional  nodal  points 
fixed  against  X-movement,  Ns  2  to  NOX ,  are 
specified  in  the  next  columns  and  on  additional 
cards  as  required. 


Explanation 


Nodal  point  number  of  the  first  nodal  point  fixed 
against  both  X-  and  Y-movement. 


Additional  nodal  points  fixed  against  both  X- 
and  Y-movement,  N  =  2  to  NOXY ,  are 
specified  in  the  next  columns  and  on  additional 
cards 

as  required. 


One  card  is  supplied  for  each  two-dimensional  element  and  interface  element;  bar 
elements  are  not  included  in  this  series  of  cards.  All  interface  elements  are  supplied  in 
sequential  order  first,  followed  by  two-dimensional  elements,  also  in  sequential  order. 
Thus,  interface  elements  must  be  numbered  from  N  =  1  to  NUMJT ,  and  two-dimen¬ 
sional  elements  from  N  =  (NUMJT  +  1)  to  NUMEL  . 

Nodal  point  numbers  must  be  specified  consecutively,  proceeding  counterclockwise 
around  the  element,  as  shown  in  Figure  A2. 

The  first  and  last  nodal  point  numbers  specified  for  interface  elements  must  have  the 
same  coordinates,  as  shown  in  Figure  A2.  Triangular  two-dimensional  elements  having 
four  different  nodal  point  numbers  may  not  be  used;  the  first  and  last  point  numbers  of 
a  triangular  element  must  be  identical. 
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Figure  A2.  Nodal  point  numbering  of  isoparametric 
QM5element  and  interface  element  (Clough  1984) 


Element  numbers  in  a  row  may  be  omitted,  in  "which  case  the  omitted  elements  will  be 
generated  by  incrementing  the  element  number  N  and  the  nodal  point  numbers  I ,  J , 
K  ,  and  L  by  one,  and  by  assigning  the  same  material  type  number  as  specified  for  the 
last  element.  The  first  and  last  elements  in  the  row  must  be  specified  (see  Figure  Al). 
If  no  elements  are  omitted,  the  element  numbering  may  be  done  in  any  order,  provided 
all  interface  elements  are  numbered  first. 


Variable 


N 


Element  number. 


IL(N,1)  Number  of  nodal  point  I . 

IL(N,2)  Number  of  nodal  point  J  . 

IL(N,3)  Number  of  nodal  point  K  . 
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Variable 


Explanation 


IL(N,4)  Number  of  nodal  point  L . 

EL(N,5)  Material  type  number. 

10.  BAR  ELEMENT  CARD 


One  card  is  supplied  for  each  bar  element  initially  in  the  mesh  or,  added  during  an  in¬ 
cremental  load  step.  Note  that  for  a  continuation  analysis  this  card  is  not  automatically 
generated.  Elements  are  numbered  sequentially  from  N  =  1  to  NMBAR  ,  with 
NSTRAN  strain  gage  bar  elements  numbered  first. 

Variable  _ Explanation _ 


N 

IB(N.l) 

IB(N,2) 

IB(N,3) 


BAR(N.l) 

BAR(N,2) 


Bar  element  number. 

Number  of  nodal  point  I . 

Number  of  nodal  point  J . 

Spring  response  type  code  - 
-  1  if  both  compression  and  tension  of  bar  allowed. 
=  2  if  only  compression  allowed. 

=  3  if  only  tension  allowed. 

cos  a 

sin  a 

The  sign  convention  of  angle  a  is  determined 
as  shown  in  Figure  A3.  Angle  a  is  measured 
counter-clockwise  from  a  line  drawn  in  the 
positive  x-direction,  originating  at  node  I  and 
connecting  node  I  to  node  J . 


BAR(N,3)  Prestress  force  in  the  bar  element.  The  force  must 

be  prescribed  at  both  node  I  and  node  J  using  a— 
load/construction  mode  9  load  step. 

BAR(N,4)  Stiffness  of  bar  element.  This  is  usually  computed 

as  AE/L  ,  but  the  mesh  length  (distance  from 
node  I  to  node  J )  need  not,  and  usually  does  not, 
correspond  to  the  actual  length.  Set  equal  to  1.0 
for  bar  elements  used  as  strain  gages. 

BAR(N,5)  Displacement  of  bar  element  at  activation.  This  al¬ 

lows  for  a  specified  magnitude  of  slack  at  the  strut 
connection;  the  bar  will  deform  BAR(N,5)  prior 
to  its  stiffness  being  activated. 

Bar  elements  can  function  as  either  anchors  or  strut  (spring)  supports.  The  required 
parameters  are  dependent  on  the  type  of  bar  element  specified. 
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If  a  strut  support  is  specified,  nodal  point  J  is  a  node  fixed  against  x-  or  y-movement, 
depending  on  the  orientation  of  the  strut  being  modeled.  For  program  storage  efficiency 
the  number  of  node  J  should  be  as  close  as  possible  to  the  number  of  node  I .  Nodal 
point  I  represents  the  point  of  connection  between  the  wall  and  the  actual  strut.  Nodal 
points  I  and  J  ,  then,  are  not  necessarily  physically  connected,  sinCe  the  element  stiff¬ 
ness  is  input  independently.  Nodal  point  J  allows  the  force  at  the  J  node  to  be  carried 
into  the  system  as  a  reaction  at  a  fixed  node.  This  is  consistent  with  the  typical  mesh 
representation  of  one  half  of  a  symmetric  excavation.  The  values  of  sin  a  and  cos  a 
are  specified  according  to  the  sign  convention  shown  in  Figure  A3.  The  values  are 


sin  a  -  - 

cos  a  •  ♦  sin  •  *  ♦ 

cos  a  -  - 


Figure  A3.  Sign  convention  and  definition  of  local  axes 
for  bar  element  (Clough  1984) 
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input  to  represent  the  line  of  action  of  the  strut  support,  and  do  not  need  to  correspond 
to  the  relative  positions  of  the  I  and  J  nodes. 

If  an  anchor  is  specified,  nodal  points  I  and  J  physically  represent  two  ends  of  the 
anchor,  and  must  be  restrained  appropriately.  The  values  of  sin  a  and  cos  a  must  cor¬ 
respond  to  the  relative  positions  of  the  I  and  J  nodes  representing  the  ends  of  the 
anchor.  Stiffness  is  computed  as  AE/L  ,  with  L  being  the  distance  between  nodes  I 
and  J ,  and  either  A  or  E  altered  to  give  the  correct  stiffness.  Stiffness  of  an  anchor 
or  a  strut  support  is  inputted  as  force  per  length  of  wall. 

For  either  element  type,  specifying  the  prestress  force  does  not  apply  the  force.  The 
concentrated  force  load/construction  mode  must  be  used  for  this  Durpose.  Thus,  bar  ele¬ 
ments  initially  in  the  mesh  cannot  carry  a  prestress  force,  since  it  is  not  applied  by  a 
gravity  turn-on  analysis. 

11.  MOMENT  COMPUTATION 

This  input  is  supplied  only  if  NOMOM  *  0 .  One  card  is  supplied  for  each  pair  of 
bar  element  strain  gages  used  in  moment  computations  (see  Equation  30). 


Variable 

MOMBAR(N.l) 


Explanation 


MOMBAR(N,  1)  Bar  element  number  for  strain  gage  No.  1. 

MOMB  AR(N,2)  Bar  element  number  for  strain  gage  No.  2. 

This  card  is  repeated  for  N  =  2  ,  NOMOM  pairs  of  bar  element  strain  gages. 

12.  SOIL  REINFORCEMENT 

This  input  is  supplied  only  if  R_FILL  *  0  .  This  series  of  cards  define  the  layer 
geometry  and  material  properties  for  the  reinforcement  material  placed  within  a  back¬ 
fill.  Each  layer  of  reinforcement  is  modeled  using  a  series  of  bar  elements.  The  total 
number  of  reinforcement  bar  elements  and  coordinates  is  determined  within  the  com¬ 
puter  code.  The  length  of  each  bar  element  corresponds  to  the  length  of  each  reinforce¬ 
ment  layer  that  is  contained  within  each  two-dimensional  soil  element.  The  number  of 
nodal  points  for  the  mesh  is  not  increased  by  the  presence  of  the  reinforcement  bar  ele¬ 
ments  because  an  embedment  technique  is  used  to  account  for  the  increased  stiffness  of 
the  backfill  as  a  result  of  the  presence  of  the  layer(s)  of  reinforcement  (see  paragraph 
29).  The  only  restriction  imposed  by  the  embedment  procedure  is  the  number  of  in¬ 
dividual  bar  elements  created  by  SOILSTRUCT.  A  bar  element  is  created  for  each 
solid  element  intersected  by  a  reinforcement  layer.  At  present,  up  to  500  bar  elements 
can  be  created. 

There  are  two  restrictions  regarding  the  layout  of  reinforcement  layers);  (1)  the  layers 
may  not  be  coincident  with  the  horizontal  sides  of  any  two-dimensional  soil  element, 
and  (2)  the  layers)  may  not  intersect  any  nodal  point  defining  a  triangular  soil  element. 
Where  these  restrictions  might  apply,  the  coordinates  of  the  reinforcement  and  the  ele¬ 
ment  side/nodes  need  only  differ  by  a  very  small  amount  ( >10r5 ). 


/ 

I 


/ 


J 


(a)  FIRST  CARD 

Variable 

Explanation 

REBARS 

Number  of  layers  of  reinforcement  material. 

BAR.STIF 

Product  of  cross-sectional  area  and  Young’s 

moduli,  (A-E),  per  unit  width. 

(b)  SECOND  CARD 

Variable 

Explanation 

MATR_FILL(N)  Material  type  number  for  two-dimensional  soil  ele¬ 
ment  containing  reinforcement. 


Additional  material  type  numbers,  N  =  2  to 
R_FILL  ,  are  specified  in  the  next  columns  and  on 
additional  cards  as  required. 


(c)  THIRD  CARD 

Variable 

Explanation 

J 

Number  of  reinforcement  layer. 

XI(J) 

X-coordinate  of  the  left-most  end  of  reinforcement 
layer  J. 

YI(J) 

Y-coordinate  of  the  left-most  end  of  reinforcement 
layer  J . 

XJ(J) 

X-coordinate  of  the  right-most  end  of  reinforce¬ 
ment  layer  J . 

YJ(J) 

Y-coordinate  of  the  right-most  end  of  reinforce¬ 
ment  layer  J . 

This  card  is  repeated  for  J  =  2  to  REBARS  layers  of  reinforcement. 

13,  CONTINUATION  OR  INITIALIZATION  CARDS 

This  input  is  supplied  only  if  INIT  =  0 ,  and  is  supplied  from  disc  storage  that  has 
been  generated  during  a  preceding  analysis.  This  option  is  provided  so  that  a  required 
sequence  of  loading  steps  can  be  stopped  at  an  intermediate  step,  then  restarted  from 
that  step  without  repeating  the  complete  analysis.  These  cards  may  also  be  used  to 
specify  values  for  particular  variable(s)  in  an  initial  analysis  without  using  the  gravity 
turn-on  procedure.  Similarly,  values  assigned  to  specific  variables  can  be  changed  if 
the  sequence  of  loading  is  stopped  prior  to  a  restart  analysis. 
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(a)  FIRST  CARD  (ELEMENT  INFORMATION)  —  FORMAT  (415) 


Column 

Variable 

Explanation 

1-5 

NUMEL 

Number  of  elements  in  the  mesh,  excluding  bar  ele¬ 
ments,  but  including  interface-elements. 

6-10 

NUMJT 

Number  of  interface  elements. 

11-15 

NUMBAR 

Number  of  bar  elements,  including  those  initially 
in  the  mesh  and  those  added  in  previous  loading 
steps  (if  a  restart  analysis). 

16-20 

NUMNP 

Number  of  nodal  points. 

(b)  SECOND  CARD  (STRESS  INFORMATION)  —  FORMAT  (4F20.6) 

One  card  is  supplied  for  each  interface  element  or  two-dimensional  element,  in  numeri- 

cal  sequence  N  =  1  to  NUMEL . 

Column 

Variable 

Explanation 

1-20 

SIG(N,1) 

Stress  in  the  x-direction  for  a  two-dimensional 
element,  ox,  normal  stress  for  an  interface  element. 

21-40 

SIG(N,2) 

Stress  in  the  y-direction  for  a  two-dimensional 
element,  oy ,  shear  stress  for  an  interface  element. 

41-60 

SIG(N,3) 

x-y  shear  stress  for  a  two-dimensional  element 
txy ,  zero  or  blank  for  interface  element. 

61-80 

SIG(N,4) 

Maximum  previous  value  of  x-y  shear  stress  for  a 
two-dimensional  element;  zero  or  blank  for  inter¬ 
face  element. 

(c)  THIRD  CARD  (NODAL 

POINT  INFORMATION)  —  FORMAT  (2D  15.9, 

1F9.2,  2X,  2D15.9,  1F9.2) 

Information  for  two  nodal  points  is  supplied  on  each  card.  Nodal  points  are  specified  in 
numerical  order,  N  =  1  to  NUMNP  . 

Column 

Variable 

Explanation 

1-15 

DISPX(N) 

Displacement  in  x-direction. 

16-30 

DISPY(N) 

Displacement  in  y-direction. 

31-39 

PP(N) 

Pore  pressure  in  head  of  water. 

41-56 

DISPX(N+1) 

Displacement  in  x-direction. 

57-71 

DISPY(N+1) 

Displacement  in  y-direction. 

72-80 

PP(N+1) 

Pore  pressure  in  head  of  water. 

Information  for  the  next  two  nodal  points  is  sup¬ 
plied  on  subsequent  lines. 
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(d)  FOURTH  CARD  (MATERIAL  TYPE  DESIGNATION)  —  FORMAT 
(1515) 

Material  type  numbers  for  15  two-dimensional  or  interface  elements  are  specified  on 
each  card.  Material  type  numbers  for  elements  in  numerical  sequence,  N  =  1  to 
NUMEL ,  are  specified. 

Column  Variable  _ Explanation _ 

1-5  IL(N,5)  Material  type  number. 

Material  type  numbers  for  the  next  14  elements  are 
supplied  in  the  next  14  five-column  fields. 

Note  that  material  type  numbers  supplied  on  these  cards  supersede  the  material  type 
numbers  specified  on  the  element  card  (section  9).  Thus  material  type  changes  can  be 
made  as  part  of  a  restart  analysis  rather  than  including  such  changes  in  a  loading  step  of 
an  analysis. 

(e)  FIFTH  CARD  (BAR  ELEMENT  INFORMATION)  —  FORMAT 
(3I5,2D10.7,  2D10.1,  Di0.5) 

This  card  is  supplied  only  if  bar  elements  are  included  (NUMB  AR  >  0).  Information 
for  bar  elements  is  specified  on  each  card.  Information  is  supplied  for  bar  elements  in 
numerical  sequence,  N  =  1  to  NUMBAR  . 


Column 

Variable 

Explanation 

1-5 

IB(N,1) 

Number  of  nodal  point  I . 

6-10 

IB(N,2) 

Number  of  nodal  point  J. 

11-15 

IB(N,3) 

Spring  response  type  code  - 
■  1  if  both  compression  and  tension  of  bar  allowed. 

2  if  only  compression  allowed, 
s  3  if  only  tension  allowed. 

16-25 

BAR(N.l) 

cos  a 

26-35 

BAR(N,2) 

sin  a 

36-45 

BAR(N,3) 

Prestress  force  in  the  bar  element. 

46-55 

BAR(N,4) 

Stiffness  of  bar  element,  (AE/L). 

56-65 

BAR(N,5) 

Displacement  of  bar  element  at  activation. 

Note  that  these  parameters,  if  changed  for  a  restart  analysis,  supersede  those  specified 
on  the  Bar  Element  Card  (section  10).  Also,  if  bar  elements  are  initially  included  in  the 
mesh,  and  this  initialization  procedure  is  used,  then  this  card  must  be  included,  duplicat¬ 
ing  the  information  specified  in  the  Bar  Element  Card. 
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(f)  SIXTH  CARD  (INTERFACE  INFORMATION)  —  FORMAT  (8D10.4) 

This  card  is  supplied  only  if  interface  elements  are  included  (NUMJT  >  0).  Information 
for  four  interface  elements  is  specified  on  each  card.  Information  is  supplied  for  inter¬ 
face  elements  in  numerical  sequence,  N  =  1  to  NUMJT . 


Column 


11-20 


Variable 
STF:  (N) 
STFN(N) 


_ Explanation _ 

Shear  stiffness  of  first  interface  element. 

Normal  stiffness  of  first  interface  element. 

Information  for  the  next  three  interface  elements  is 
supplied  in  the  next  six  10-column  fields. 


Note  that  the  value  of  the  shear  stiffness,  if  changed  for  a  restart  analysis,  supersedes 
the  value  specified  on  the  Interface  Property  Card  (section  6).  Thus  the  interface  stiff¬ 
ness  can  be  changed  as  part  of  a  restart  analysis. 

(g)  SET  OF  5  CARDS  (REINFORCEMENT  BAR  INFORMATION) 

This  set  of  5  cards  is  supplied  only  if  reinforcement  bar  elements  are  included 
(R_FILL  *  O) . 

Card  1  of  5  —  FORMAT  (15) 


Column 


Column 


Variable 

REBAR.NUM 


_ Explanation 

I 

Number  of  reinforcement  bar  elements. 


Card  2  of  5  —  FORMAT  (1515) 


Variable 

EMBED(N.l) 


Explanation 


EMBED(N.l)  Two-dimensional  soil  element  number  in  which 

the  reinforcement  bar  number  N  is  contained. 

! 

Information  for  N  =  1  to  REBAR_NUM  reinfor¬ 
cement  bar  elements  are  to  be  provided.  The  next 
14  EMBED(N,1)  values  are  supplied  in  the  next 
14  columns. 

Card  3  of  5  —  FORMAT  (1D15.9,  1025.19,  1D15.9,  1D25.19) 


Column 


16-40 


41-55 

56-80 


Variable 

R_COORD(N,  1,1) 
R_COORD(N,2, 1 ) 

R_COORD(N,  1 ,2) 
R_COORD(N,2,2) 


Explanation 


X  coordinate  for  reinforcement  bar  element  node  I. 

Local  X  coordinate  for  reinforcement  bar  element 
node  T. 

X  coordinate  for  reinforcement  bar  element  node  J. 

Local  X  coordinate  for  reinforcement  bar  element 
node  J. 


There  are  REBAR_NUM  cards  provided. 
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Card  4  of  5  —  FORMAT  (8D10.4) 


Column 

Variable 

Explanation 

1-10 

STIF_TEMP(N) 

Product  of  cross-sectional  area  and  Young’s 
moduli,  (A-E) .  per  unit  width  for  reinforceme 
bar  element  av . 

Information  for  the  next  seven  reinforcement 
elements  is  supplied  in  the  next  seven  fields 

There  are  REBAR_NUM  values  provided. 

Card  5  of  5  —  FORMAT  (7D 1 1 .5) 

Column 

Variable 

Explanation 

1-11 

R_FORCE(N,2) 

Reinforcement  bar  element  strain. 

Information  for  the  next  six  reinforcement  ele¬ 
ments  is  supplied  in  the  next  six  fields.  There  are 
REBARJNUM  values  provided. 


14.  LOADING  STEP  CARDS 

These  cards  are  required  only  if  the  number  of  load  steps  is  greater  than  zero,  NC  *  0 . 
Up  to  three  load/construction  modes  may  be  used  per  load  step,  as  described  in  section  3. 
Cards  are  assembled  in  the  order  specified  on  the  Loading  Information  Cards  (section  3): 
cards  for  the  first  load/construction  mode  specified  for  the  first  loading  step,  and  so  on, 
to  the  cards  for  the  last  load/construction  mode  specified  for  the  last  construction  step. 

For  a  given  loading  step,  the  lowest  numbered  load/construction  mode  is  processed 
first,  but  the  analysis  of  the  loading  step  is  made  for  the  combined  effect  of  all  load/ 
construction  modes  specified  for  that  loading  step.  Care  must  therefore  be  exercised 
in  specifying  some  load/construction  modes,  such  as  material  type  changes  or  con¬ 
centrated  forces  representing  prestress  forces,  in  the  same  loading  step  with  other 
load/construction  modes. 

(a)  EXCAVATION 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(n,2) ,  or  KCS(N,3)  =  1 .  Input  is 
handled  by  subroutine  EXCAV  .  Free-excavated  elements  and  common-excavated 
elements  are  input  separately.  A  free-excavated  element  is  an  element  specified  t'  be 
excavated  that  has  no  boundary  in  common  with  an  element  not  specified  to  be  ex 
cavated  in  the  loading  step.  A  common-excavated  element,  therefore,  has  at  least  te 
boundary  in  common  with  an  unexcavated  element  (Figure  A4). 

Interface  elements  can  only  be  included  as  free-excavated  elements,  even  if  they  ha  a 
boundary  in  common  with  an  unexcavated  element.  Interface  etetuwut s  cannot  be  us  i 
as  interpolation  elements.  Free-excavated  elements  (other  than  interpolation  elements) 
can  be  used  as  interpolation  elements,  though  common-excavated  elements  are  more 
commonly  used. 
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ELEMENTS  1  AND  2  ARE  FREE  EXCAVATED  ELEMENTS. 
ELEMENTS  3.  6.  7.  AND  8  ARE  COMMON  EXCAVATED 
ELEMENTS. 

NODES  13.  14.  15.  16.  22.  AND  28  ARE  LOADED 
BY  EXCAVATION  FORCES. 

ELEMENTS  6.  7.  AND  8  SHOULD  BE  INPUT 
SEQUENTIALLY  FOR  OPTIMUM  EFFICIENCY. 


Figure  A4.  Example  excavation  load  step  defining  free-  and 
common-excavated  elements  in  relation  to  the  excavation 
boundary  (Clough  1984) 

If  possible,  adjacent  common-excavated  elements  should  be  input  sequentially  as  this 
procedure  avoids  repetitive  computation;  nodal  loads  for  a  nodal  point  common  to  the 
two  sequential  elements  will  only  be  calculated  once. 

(1)  FIRST  CARD  (CONTROL  DATA)  —  FORMAT  (215) 

Column  Variable  _ Explanation _ 

1-5  NFXEL  Number  of  free-excavated  elements. 

6-10  NXELCB  Number  of  common-excavated  elements. 

(2)  SECOND  CARD  (FREE-EXCAVATED  ELEMENT  DATA)  — 
FORMAT  (1615) 

Element  numbers  of  16  free-excavated  elements  can  be  supplied  on  one  card.  A  maxi¬ 
mum  of  150  can  be  specified  in  one  loading  step.  Element  numbers  of  all  free  ex¬ 
cavated  elements,  N  =  1  to  NFXEL ,  are  to  be  specified. 


/ 

'J- 
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Column  Variable  _  Explanation _ 

1-5  LNXEL(N)  Element  number  of  first  free-excavated  element. 

Information  for  the  next  15  free-excavated  elements 
is  supplied  in  the  next  15  five-oolumn  fields. 

(3)  THIRD  CARD  (COMMON-EXCAVATED  ELEMENT  DATA)  — 
FORMAT  (815) 

One  card  is  supplied  for  each  common-excavated  element,  N  =  1  to  NXELCB  .  A 
maximum  of  150  common-excavated  elements  can  be  specified  in  one  loading  step. 
Loading  codes  include: 

0  •  the  node  is  not  loaded  by  excavation  forces,  and  is  not  common  to  both  an 
excavated  and  an  unexcavated  element. 

1  -  the  node  is  loaded  by  excavation  forces  and  is  common  to  both  an  excavated 
and  an  unexcavated  element. 

Notes  I ,  J ,  K ,  and  L  refer  to  the  same  nodes  I ,  J  ,  K  ,  and  L  specified  on  the 
element  card  (section  9). 

Interpolation  elements  should  be  specified  in  a  crisscross  fashion,  as  shown  in  Figure  A5. 
Further,  the  x-  or  y-coordinates  of  diagonal  elements  must  not  be  the  same.  If  these 
rules  are  not  adhered  to,  the  interpolation  routine  will  detect  a  singularity,  and  process¬ 
ing  will  stop. 


Column 

Variable 

Explanation 

1-5 

LUL(N.l) 

Element  number  of  the  first  common-excavated 
element.  This  element  is  also  the  first  interpola¬ 
tion  element. 

6-10 

LUL(N,2) 

Element  number  of  second  interpolation  element. 

11-15 

LUL(N,3) 

Element  number  of  third  interpolation  element. 

16-20 

LUL(N,4) 

Element  number  of  fourth  interpolation  element. 

21-25 

LUL(N,5) 

Loading  code  for  node  I . 

26-30 

LUL(N,6) 

% 

Loading  code  for  node  J . 

31-35 

LUL(N,7) 

Loading  code  for  node  K  . 

36-40 

LUL(N,8) 

Loading  code  for  node  L . 

Loading  codes  are  specified  for  the  nodes  of  the 
element  to  be  excavated  specified  in  columns  1-5. 

(b)  FILL  OR  CONCRETE  PLACEMENT 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =  2  .  Input  is 
handled  by  subroutine  BUILD  .  The  types  of  elements  that  can  be  placed  include 
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structural,  soil,  and  interface  elements.  Subroutine  SUBSTP  cannot  be  used  in  con¬ 
junction  with  subroutine  BUILD . 

At  placement,  the  fill  element(s)  is  assigned  a  low  modulus  and  the  surface  displace¬ 
ments  are  set  equal  to  zero.  Stresses  assigned  to  the  newly  placed  fill  are  based  on 
ox  *  Kjj  oy  or  a%  =  oy  (see  KEYH20  on  card  grour  .  wheic  ^  -Huai  to  the 

product  of  the  unit  weight  of  the  fill  element  times  the  depth  below  the  surface  to  the 
center  of  the  element.  If  a  fill  element  contains  a  reinforcement  layer(s),  then  the  bar 
element(s)  representing  the  layerfs)  will  be  added  to  the  newly  placed  fill  elements  at 
the  end  of  the  load  step.  Generally,  all  fill  elements  associated  with  a  particular  layer 
should  be  added  in  one  load  step. 

(1)  FIRST  CARD  (CONTROL  DATA)  —  FORMAT  (215,  D10.2, 15) 


Column 

1-5 


Variable 

NLEL 


Explanation 


Total  number  of  elements  to  be  placed,  including 
interface  and  structure  elements. 

6-10  NONP  |  Number  of  nodal  points  within  layerfs)  to  be  as¬ 

signed  zero.  This  includes  all  nodal  points  for 
existing  elements  to  be  placed  except  points  in 
j  common  with  an  existing  element. 

1 1-20  HTB  New  y-coordinate  at  the  top  of  backfill. 

I 

21-25  NPANEL  |  Number  of  segment  end  points  used  to  define  the 

y-coordinates  for  the  top  of  the  backfill.  NPANEL 
must  be  less  than  or  equal  to  30.  If  NPANEL  «=  0 , 
elevation  of  backfill  =  HTB  value. 

(2)  SECOND  CARD  (ELEMENT  NUMBER)  —  FORMAT  (1615) 

Element  numbers  of  16  placell  elements  can  be  supplied  on  the  card.  A  maximum  of 
250  can  be  specified  in  one  loading  step.  Element  numbers  of  all  placed  elements, 

N  ■  1  to  NLEL  ,  are  to  be  supplied. 

Column  Variable  _ Explanation _ 


1-5  LEL(N,1)  Element  number  of  first  element  to  be  placed. 

6-10  LEL(N,2)  Material  type  number. 

Information  for  the  next  seven  place  element  num¬ 
bers  and  their  corresponding  material  types  are 
supplied  in  the  next  14  five-column  fields. 

(3)  THIRD  CARD  (NODE  NUMBERS)  —  FORMAT  (1615) 

Nodal  point  numbers  of  16  nodes  to  be  assigned  zero  displacement  can  be  supplied  on 
one  card.  A  maximum  of  250  can  be  specified  in  one  loading  step.  Nodal  point  num¬ 
bers,  N  *  1  to  NONP ,  are  to  be  specified  in  sequential  order. 
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Column 

Variable 

Explanation 

1-5 

NP(N) 

Nodal  point  number  of  first  node  assigned  zero  dis¬ 
placement. 

Information  for  the  next  15  nodal  points  is  sup¬ 
plied  in  the  next  15  five-column  fields. 

(4)  FOURTH  CARD 

(TOP  OF  BACKFILL)  —  FORMAT  (2D10.2) 

Column 

Variable 

Explanation 

1-10 

XPANEL(N) 

X-coordinate  of  ground  surface  segment  end  point. 

11-20 

YPANEL(N) 

Elevation  of  ground  surface  at  the  segment  end 
point. 

Information  for  the  remaining  segment  end  points 
is  supplied  on  subsequent  cards. 

(c)  SEEPAGE 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =  3  .  Input  is  hand¬ 
led  by  subroutine  SEEP .  Seepage  loads  are  determined  from  change  in  pore  pressure 
specified  as  DP(N)  on  nodal  point  cards,  or  from  the  specified  phreatic  level  changes. 

(1)  FIRST  CARD  (CONTROL  DATA)  —  FORMAT  (15) 

Column 

Variable 

Explanation 

1-5 

NCODE 

Option  code  specifying  how  seepage  loading  data 
is  to  be  input  - 

»  0  if  specified  as  DP(N)  on  Nodal  Point  Card 
(section  7). 

=  1  if  the  seepage  loading  is  to  be  calculated 
using  the  new  phreatic  surface  input  on  the 
following  cards. 

— 

(2)  SECOND  CARD  (NUMBER  OF  PHREATIC  SEGMENTS)  — 
FORMAT  (15) 

This  card  is  required  only  if  NCODE  =  1  . 

Column 

Variable 

Explanation 

1-5 

NWAT 

Number  of  phreatic  surface  segment  end  points 
used  to  specify  the  new  phreatic  surface.  NWAT 
must  be  greater  than  or  equal  to  2.  The  number  of 
phreatic  surface  segments  is  equal  to  NWAT . 

The  maximum  value  of  NWAT  is  30. 

(3)  THIRD  CARD  (PHREATIC  LEVEL  DATA)  —  FORMAT  (6D10.2) 
This  card  is  required  only  if  NCODE  =  1  . 
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The  end  points  of  the  phreatic  surface  segments  delineating  the  new  and  old  phreatic 
surfaces  are  specified  as  x-coordinates  and  must  be  the  same  as  the  x-coordinate  of  a 
nodal  point  (Figure  A6).  Both  the  present  and  new  phreatic  surfaces  are  assumed  to  be 
linear  between  the  bounding  x-coordinates.  The  left-hand  side  of  the  mesh  is  always 
the  first  bounding  x-coordinate  specifi'd.  A  bounding  x-coordinate  on  the  old  phreatic 
surface  will  require,  usually,  specification  of  the  same  x-coordinate  on  uie  new  phreatic 
surface. 

Two  end  points  (x-coordinate),  with  associated  new  and  old  phreatic  levels  (y-coor- 
dinates),  are  supplied  on  each  card.  All  end  points,  N  -  1  to  NWAT ,  must  be 
specified. 

Column  Variable  _ Explanation _ 

I- 10  XW(N)  X-coordinate  bounding  the  levels  PREL(N)  and 

FUEL(N)  on  the  right-hand  side.  (Must)  be  the 
same  as  the  x-coordinate  of  a  nodal  point. 

II- 20  PREL(N)  Present  level  (y-coordinate)  of  the  phreatic  sur¬ 

face  at  XW(N) . 


21-30  FUEL(N)  New  level  (y-coordinate)  of  the  phreatic  surface  at 

XW(N) . 

Information  for  next  end  point  is  supplied  in  the 
next  three  10-column  fields. 
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Figure  A6 .  Example  illustrating  phreatic  level  data  required 
for  the  seepage  loading! construction  case  ( Clough  1984) 


A30 


(d)  DELETION  OF  BAR  ELEMENTS 

These  cards  are  supplied  only  if  KCS(N,1)  ,  KCS(N,2) ,  or  KCS(N,3)  =  4  .  Input  is 
handled  by  the  main  program. 

The  deleted  bar  elements  remain  in  the  mesh  but  with  zero  stiffness.  The  force  the  bar 
element  carried  is  applied  to  the  free  node  or  nodes  at  its  ends. 

This  loading/construction  mode  cannot  be  specified  in  the  same  loading  step  as  fill  or 
concrete  placement. 

(1)  FIRST  CARD  (CONTROL  DATA) 

Column  Variable  _ Explanation _ 

1-5  NCARDS  Number  of  deleted  bar  elements.  There  is  no  limit 

other  than  the  number  of  bar  elements  presently  in 
the  mesh. 

(2)  SECOND  CARD  (ELEMENT  NUMBERS) 

The  element  numbers  of  16  deleted  bar  elements  can  be  specified  on  one  card.  A  total 
of  N  =  1  to  NCARDS  cards  must  be  supplied. 


Column 

Variable 

Explanation 

1-5 

N 

Element  number  of  bar  element  to  be  deleted. 

Element  numbers  for  the  next  15  elements  are  sup¬ 
plied  in  the  next  15  five-column  fields. 

(e)  ADDITION  OF  BAR  ELEMENTS 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =  5  .  Input  is 
handled  by  the  main  program. 

Information  on  the  second  card  is  the  same  as  that  explained  for  the  Bar  Element  Card 
(Section  10).  The  added  bar  elements  are  numbered  sequentially  from  NUMBAR  +  1 , 
where  NUMBAR  is  the  number  of  bar  elements  in  the  mesh  before  the  present  loading 
step. 

(1)  FIRST  CARD  (CONTROL  DATA) 

Column  Variable  Explanation _ 

1-5  NCARDS  Number  of  bar  elements  to  be  added.  Any  number 

can  be  added  in  a  loading  step;  however,  the  maxi¬ 
mum  number  of  bar  elements  that  can  be  in  the 
mesh  (including  inactive  or  deleted  elements)  is  15. 

(2)  SECOND  CARD  (ADDED  BAR  ELEMENT  DATA) 

Variable  _ _ Explanation _ 

N  Element  number  of  added  bar  element. 

IB(N,1)  Number  of  nodal  point  I . 
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Variable 

Explanation 

IB(N,2) 

Number  of  nodal  point  J . 

IB(N,3) 

Spring  response  type  code. 

BAR(N.l) 

cos  a 

BAR(N,2) 

sin  a 

BAR(N,3) 

Prestress  force  in  the  bar  element. 

BAR(N,4) 

Stiffness  of  bar  element. 

BAR(N,5) 

Displacement  of  bar  element  at  activation. 

Note  that  in  a  restart  analysis  the  added  bar  element(s)  information  must  be  included  in 
the  bar  element  connectivity  cards  (Section  10)  and  the  value  for  NUMBAR  in  the 
data  control  card  (Section  2a)  must  be  updated. 

(f)  BOUNDARY  PRESSURE  LOADING 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =»  6 .  Input  is 
handled  by  subroutine  SURFLD  .  Linear  pressure  distributions  are  assumed,  based  on 
the  pressures  specified  for  the  nodal  points. 

(1)  FIRST  CARD  (CONTROL  DATA)  —  FORMAT  (15) 

Column  Variable  _ Explanation _ 

I- 5  NLDS  Number  of  loaded  boundaries. 

There  is  no  limit  to  the  number  of  loaded  boun¬ 
daries  that  can  be  specified. 

6-10  NOLDSX  Number  of  loaded  boundaries  for  which  the 

horizontal  components  of  the  nodal  points  load 
vectors  are  set  equal  to  zero. 

II- 15  NOLDSY  Number  of  loaded  boundaries  for  which  the  verti¬ 

cal  components  of  the  nodal  points  load  vectors 
are  set  equal  to  zero. 

(2)  SECOND  CARD  (LOADED  BOUNDARY  DATA)  —  FORMAT  (215, 
4D10.2) 

Information  for  one  loaded  boundary  is  specified  on  each  card.  Nodes  I  and  T  are 
specified  counterclockwise  on  an  element  (Figure  A7).  A  normal  cpmpressh  Taction 
(pressure)  is  positive.  A  tangential  traction  (shear  stress)  is  positive!' when  directed 
clockwise  (node  J  to  node  I )  as  shown  in  Figure  A7.  A  total  of  N  =  1  to  ('  LDS 
+  NOLDSX  +  NOLDSY)  loaded  boundaries  must  be  specified.  NLDS  cards  are 
specified  first,  followed  by  NOLDSX  cards  and  then  NOLDSY  cards. 


Nodal  point  number  of  the  first  node  of  the  loaded 
boundary. 

Nodal  point  number  of  the  second  node  of  the 
loaded  boundary. 
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Column 


Variable 


Explanation 


11-20 

21-30 

31-40 

41-50 


Value  of  the  normal  pressure  acting  at  node  I 
Value  of  the  normal  pressure  at  node  J  . 
Value  of  the  shear  stress  acting’  at  node  I . 
Value  of  the  shear  stress  acting  at  node  J  . 


(g)  TEMPERATURE  LOADING 

No  cards  are  required  for  this  load/construction  mode  (Temperature  Card,  section  7). 
If  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  «  7  .  Then  the  values  of  DP(N)  are  acknow¬ 
ledged  by  the  main  program  and  processed  as  temperature  changes.  Note  that  DP(N) 
can  also  be  used  to  input  phreatic  level  changes  for  the  seepage  loading/construction 
mode.  Thus,  if  seepage  is  specified  as  being  input  through  values  of  DP(N)  ,  seepage 
and  temperature  loading  cannot  be  included  in  the  same  loading  step.  Generally, 
temperature  loading  requires  a  restart  analysis,  with  the  DP(N)  values  being  changed 
to  reflect  the  temperature  changes  prior  to  the  analysis. 

The  temperature  scale  used  (°C  or  °F)  must  correspond  to  the  coefficient  of  thermal  ex¬ 
pansion  designated  on  the  Material  Property  Cards  (Section  5).  Temperature  changes 
are  typically  designated  for  structural  materials  only. 

(h)  CONCENTRATED  FORCE  OR  DISPLACEMENT  LOADING 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =  8  or  9  . 

Input  is  handled  by  the  main  program.  Refer  to  the  Loading  Information  Card  (section  3) 
for  instructions  on  using  loading/construction  modes  8  and  9. 

(1)  FIRST  CARD  (CONTROL  DATA) 

Variable  _ Explanation _ 

NUMNDE  Number  of  loaded  or  displaced  nodes.  There  is  no 

limit  to  the  number  of  loaded  or  displaced  nodes 
that  can  be  specified. 

(2)  SECOND  CARD  (LOAD  OR  DISPLACEMENT  DATA) 

Information  for  two  loaded  or  displaced  nodal  points  is  supplied  on  each  card.  A 
total  of  N  =  1  to  NUMNDE  nodes  must  be  specified.  Sign  convention  is  positive  to 
the  right  (positive  x-direction)  and  positive  up  (positive  y-direction).  Nodal  points 
specified  as  being  loaded  or  displaced  do  not  have  to  be  in  numerical  order. 


Variable 


_ Explanation _ 

Node  number  of  the  Erst  loaded  or  displaced  node. 

Component  of  force  or  displacement  in  the  x-direc¬ 
tion  at  node  I . 

Component  of  force  or  displacement  in  the  y-direc- 
tionatnode  I. 
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Variable 


Explanation _ 

J  Node  number  of  the  second  loaded  or  displaced 

node. 

X2  Component  of  force  or  displacement  in  the  x-direc- 

tion  at  node  J . 

Y2  Component  of  force  or  displacement  in  the  y-direc- 

tion  at  node  3  . 

Displaced  nodes  are  to  be  included  in  Boundary  Condition  Cards  (Section  8). 

(i)  ELEMENT  MATERIAL  TYPE  CHANGE 

These  cards  are  supplied  only  if  KCS(N,1) ,  KCS(N,2) ,  or  KCS(N,3)  =  10 .  Input  is 
handled  by  the  main  program. 

The  material  type  of  the  specified  element  is  changed  before  the  analysis  of  the  loading 
step  which  specifies  the  change.  The  material  type  change  includes  modifying  the 
values  of  modulus  E  and  Poisson's  ratio  GUI  and  zeroing  the  stresses  SIG(N,1) . 
Thus,  if  a  material  type  change  is  specified  in  conjunction  with  boundary  loading  (in 
the  same  loading  step)  the  elements  whose  material  type  is  changed  will  respond  to  the 
loading  with  new  material  properties. 

As  included,  this  load/construction  mode  is  intended  to  physically  represent  the  grout¬ 
ing  of  an  anchor.  At  a  given  step  in  the  analysis,  the  material  types  of  soil  elements  can 
be  changed  to  represent  the  assumed  linear  elastic  grout  zone.  If  there  is  a  need  to 
change  material  types  for  any  other  reason,  this  can  be  done  by  stopping  the  analysis 
after  the  appropriate  load  step,  modifying  the  material  types  on  the  material  type  desig¬ 
nation  card  (Section  13d),  then  restarting  the  analysis. 

(1)  FIRST  CARD  (CONTROL  CARD) 

The  maximum  number  of  elements  whose  material  type  number  can  be  changed  in  a 
load  step  is  130.  The  excavation  and  material  type  change  loading/construction  modes 
cannot  be  specified  in  the  same  loading  step  since  the  same  variable,  LUL(N,I) ,  is 
used  to  input  data  for  both. 

Variable  _ Explanation _ 

NELCH  Number  of  elements  whose  material  type  number 

is  being  changed. 

(2)  SECOND  CARD  (ELEMENT  DATA) 

A  total  of  N  =  1  to  NELCH  element  numbers  and  new  material  type  numbers  must  be 
specified. 

Variable  _ Explanation _ 

LUL(N,7)  Element  number  of  first  element  with  a  specified 

new  material  type  number. 
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Variable 


Explanation 


LUL(N,8)  New  material  type  number  of  the  specified  ele¬ 

ment. 

Information  for  the  subsequent' elements  is  sup¬ 
plied  in  the  next  columns  anu  Uu  auuiuonal  cards 
as  required. 

Output  Data  Sequence 

The  input  data  read  in  the  main  program  is  tinted  out  prior  to  the  analysis  of  the 
first  load  step.  Nodal  point  data,  element  data,  and  initial  stresses  are  printed  out  only 
if  specified.  Loading/construction  mode  information  for  a  load  step  is  printed  out  prior 
to  the  analysis  of  the  load  step.  Analysis  results,  excluding  nodal  displacement,  can  be 
specified  tc  be  printed  out  after  each  iteration  or  after  the  final  iteration  of  each  substep 
of  a  loading  step.  Continuation  data  for  a  restart  analysis  can  be  stored  on  disks  (in  a 
card  image  format),  or  punched  on  cards,  after  the  final  iteration  of  the  analysis  of  the 
last  substep  of  the  last  loading  step.  Nodal  displacements  can  be  stored  on  disks  (in  a 
card  image  format),  or  punched  on  cards,  after  the  final  iteration  of  the  analysis  of  each 
substep  of  a  loading  step. 

Results  for  all  element  types  for  all  iterations,  with  the  exception  of  the  final  itera¬ 
tion,  are  based  on  one-half  the  incremental  stress  being  used  to  update  the  previous 
stress.  Thus,  they  do  not  reflect  the  full  effect  of  the  applied  loads.  These  values  of 
stress  are  used  to  update  modulus  and  stiffness  values  while  iterating. 

Different  results,  of  course,  are  printed  out  for  the  different  types  of  finite  elements 
used.  All  output  pages  are  identified  by  the  analysis  description.  Headings  designating 
the  type  of  results  (element-type),  the  load  step  number,  the  substep  number,  and  the 
iteration  number  are  included  following  the  analysis  description. 
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1.  DISPLACEMENTS 


Displacements  for  all  nodes  are  printed  out  after  the  last  iteration  of  each  substep  of 
a  loading  step. 


Column 

Heading 

1-10 

NODAL 

POINT 

11-18 

X 

19-26 

Y 

27-40 

TOTAL 

UX 

41-54 

TOTAL 

UY 

55-68 

INCREM 

UX 

69-82 

INCREM 

UY 

83-93 

PORE 

PRESS 

94-100 

NODAL 

POINT 

_ Explanation _ 

Nodal  point  number. 

X-coordinate  of  the  node. 

Y-coordinate  of  the  node. 

Total  displacement  of  the  node  in  the  X -direction. 

Total  displacement  of  the  node  in  the  Y-direction. 

Incremental  displacement  of  the  node  in  the 
X-direction. 

Incremental  displacement  of  the  node  in  the 
Y-direction. 

Pore  pressure  acting  on  the  node  in  the  head  of 
water. 

Nodal  point  number. 


2.  RESULTS  FOR  TWO-DIMENSIONAL  ELEMENTS 


Results  for  two-dimensional  elements  are  printed  as  requested  using  the  input  data 
described  in  section  2a.  Elements  numbered  from  N  =  (NUMJT  +  1)  to  NUMSOL 
are  included. 

(a)  FIRST  CARD:  STRESSES  &  NEW  MODULUS  VALUES 


Column 

Heading 

Explanation 

2-5 

ELE 

NO 

Element  number. 

7-14 

X 

X-coordinate  of  element  center. 

16 

E  or  T 

Letter  E  indicating  effective  stresses  output  if  a 
drained  material  or  T  indicating  total  stress  out¬ 
put  if  an  undrained  material. 

18-28 

SIGMA 

X 

Normal  stress  acting  in  the  X-direction. 

30-40 

TAU 

XY 

Shear  stress  acting  in  the  X-  and 

Y-directions. 

42-52 

SIGMA 

1 

Major  principal  stress. 
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Column 

Heading 

Explanation 

54-60 

THETA 

Angle  of  the  major  principal  plane,  in  degrees  mea¬ 
sured  counterclockwise  from  the  positive  x-axis. 

62-71 

NEWE 

New  value  for  Young’s  moduli. 

73-82 

NEW  B 

New  value  for  bulk  moduli. 

83-87' 

NEW  NU 

New  value  for  Poisson’s  ratio. 

88-94 

PORE 

PRESS 

Pore  pressure. 

(b) 

SECOND  CARD: 

STRESS  VALUES  AND  OLD  MODULUS  VALUES 

Column 

Heading 

Explanation 

7-14 

Y 

Y-coordinate  of  element  center. 

18-28 

SIGMA 

Y 

Normal  stress  acting  in  the  Y-direction. 

30-40 

TAU  MAX 
PREVIOUS 

Maximum  previous  shear  stress. 

42-52 

SIGMA 

3 

Minor  principal  stress. 

54-60 

STRESS 

RATIO 

Stress  level  or  ratio. 

62-7 1 

OLDE 

Old  value  for  Young’s  moduli. 

73-82 

OLD  BULK 

Old  value  for  bulk  moduli. 

83-87 

OLD  NU 

Old  value  for  Poisson’s  ratio. 

3.  RESULTS  FOR  INTERFACE  ELEMENTS 

Results  for  interface  elements  are  printed  out  following  two-dimensional  elements,  as 
specified  by  ITRD  .  Elements  numbered  from  N  =  1  to  NUMJT  are  included. 

Column 

Heading 

Explanation 

1-12 

ELEM  NO 

Element  number. 

13-20 

X 

X-coordinate  of  element  center. 

21-28 

Y 

Y-coordinate  of  element  center. 

29-44 

NORMAL 

STRESS 

Normal  stress  acting  on  the  element. 

45-60 

SHEAR 

STRESS 

Shear  stress  acting  on  the  element. 

61-76 

NORMAL 

STIFF 

Normal  stiffness. 
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Column 

Heading 

Explanation 

77-92 

SHEAR 

STIFF 

Shear  stiffness. 

93-10G 

SL 

Stress  level. 

101-112 

PORE 

PRES 

Pore  pressure  acting  on  the  element. 

113-122 

ELEM  NO 

Element  number. 

4.  RESULTS  FOR  BAR  ELEMENTS 

Results  for  bar  elements  are  printed  following  the  results  for  the  interface  elements. 
Bar  elements  numbered  from  N  =  1  to  NUMBAR  are  included. 

Column 

Heading 

Explanation 

1-13 

BAR  ELEM 

Bar  element  number. 

14-19 

I 

Nodal  point  number  of  the  I  node  of  the  element. 

20-25 

J 

Nodal  point  number  of  the  J  node  of  the  element. 

26-33 

TYPE 

Bar  element  response  type  code  (see  section  10  of 
the  input  listing). 

34-47 

COMPR 

Force  in  the  bar  element  (positive  if  compressive). 

48-61 

COMPRESSION 

Change  in  length  of  the  bar  element  (positive  if 
compressive). 

62-75 

STIFFNESS 

Stiffness  of  the  bar  element. 

76-85 

COSA 

cos  a. 

86-95 

SINA 

sin  a. 

96-105 

BAR  ELEM 

Bar  element  number. 

5.  RESULTS  FOR  BAR  ELEMENTS  USED  AS  STRAIN  GAGES 

The  results  for  bai;  elements  used  as  strain  gages  are  printed  following  the  results  for 
the  bar  elements.  The  bar  elements  numbered  from  N  =  1  to  NSTRAN  are  included. 


Column 

Heading 

Explanation 

1-7 

BAR  ELEM 

Element  number  for  bar  used  as  strain  gage. 

8-20 

STRAIN 

BAR  STF 

LESS 

Strain  in  bar  element. 

21-35 

Stress  in  steel  at  location  of  bar  element  used  as  a 
strain  gage. 

36-51 

CONCRETE 

STRESS 

Stress  in  concrete  at  location  of  bar 
element  used  as  a  strain  gage. 

A39 


/ 


6.  RESULTS  OF  MOMENT  COMPUTATIONS 


The  results  for  each  pair  of  bar  element  strain  gages  used  in  the  moment  computa¬ 
tions  are  printed  following  the  results  for  the  bar  elements  used  as  strain  gages. 


Column 

1-10 

11-20 

21-30 

31-45 

46-60 

61-70 

71-80 

8J-95 

96-110 

111-125 


Heading 
PAIR  No. 
X-CENTER 
Y-CENTER 
GAGE  No.  1 
GAGE  No.  2 


AXIAL  FORCE 


MOMENT 


^Yotenation 


Number  of  bar  element  strain  gage  pair. 

X-coordinate  of  center  of  member. 

Y-coordinate  of  center  of  member. 

Bar  element  number  for  strain  gage  No.  1. 

Bar  element  number  for  strain  gage  No.  2. 

Young’s  modulus  of  equivalent  member. 

Width  of  equivalent  member. 

Moment  of  inertia  of  equivalent  member. 

Resultant  force  acting  normal  to  cross-section  of 
equivalent  member. 

Moment  at  equivalent  member. 


7.  RESULTS  FOR  REINFORCEMENT  ELEMENTS 


The  results  for  bar  elements  used  as  soil  reinforcement  elements  are  printed  follow- 


ing  the  moment  computations. 

Column  Heading 

Explanation 

1-4 

BAR  No. 

Reinforcement  bar  element  number. 

5-13 

ELM  No. 

Soil  element  number  in  which  reinforcement  bar 
element  is  embedded. 

14-21 

X_i 

X-coordinate  of  bar  node  - 1 . 

• 

22-34 

YJ 

Y-coordinate  of  bar  node  I . 

35-47 

XJ 

X-coordinate  of  bar  node  J  . 

48-60 

YJ 

Y-coordinate  of  bar  node  J . 

61-74 

FORCE  per 

WIDTH 

Force  in  bar  per  unit  width. 

75-95 

AREA- 

MODULUS 

Product  of  cross-sectional  area  and  Young’s 
modulus. 

96-111 

STRAIN 

Strain  in  bar. 
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APPENDIX  B:  SEQUENCE  OF  OPERATIONS 

1.  The  program  SOILSTRUCT  uses  the  direct  stiffness  method,  {F}  =  [K] { u }  ,  to 
solve  for  incremental  nodal  displacements  { u)  resulting  from  incremental  loads  ap¬ 
plied  to  the  nodal  points  {F}  .  The  local  element  stiffness  matrices  are  first  formulated, 
then  assembled  into  the  global  stiffness  matrix.  Equivalent  nodal  loads  due  to  construc¬ 
tion  or  applied  loadings  are  assembled  in  the  incremental  load  vector,  {F}  .  The  com¬ 
puted  incremental  displacements  are  then  used  to  compute  the  incremental  change  in 
stress  acting  at  the  center  of  the  elements.  The  values  of  total  stress  are  updated  by  the 
computed  incremental  changes  in  stress.  The  total  stresses  are  then  used  to  revise  the 
elastic  moduli  used  in  the  formulation  of  the  element  stiffnesses.  These  procedures  are 
repeated  for  each  iteration,  and  in  turn  for  each  substep  of  each  load  case  during  the 
analysis.  The  determination  of  initial  gravity  stresses  is  accomplished  in  a  similar  fash¬ 
ion  and  can  be  viewed  as  an  initial  load  step,  with  the  nodal  loads  equal  to  the  body 
forces. 

2.  The  following  is  a  listing  of  the  names  of  each  of  the  26  subroutines  comprising 
the  program  SOILSTRUCT  and  a  brief  description  of  their  purpose: 

•  Main  Program.  The  main  program  serves  to  control  the  execution  of 
SOILSTRUCT.  It  calls  subroutines,  prints  input  data,  load  case  information,  ma¬ 
terial  properties,  node  and  element  data,  and  boundary  conditions.  The  input 
data  for  excavation,  seepage,  embankment  construction,  and  boundary  loadings 
are  printed  in  their  respective  subroutines.  Calculated  equivalent  nodal  loads 
due  to  installation  or  deletion  of  bar  elements,  and  concentrated  forces  or  dis¬ 
placements,  are  added  in  the  main  program. 

•  DETNA.  Subroutine  DETNA  calculates  the  number  of  degrees  of  freedom,  de¬ 
termines  the  locations  of  the  diagonal  terms  of  the  global  stiffness  matrix  in  the 
vector  SN ,  and  computes  the  required  size  of  SN  . 

•  INITAL.  Subroutine  INITAL  calculates  and  prints  initial  stresses  for  a  gravity 
tum-on  analysis.  This  is  done  by  sequential  calls  to  STRSTF,  OPTSOL,  and 
STRESS.  If  a  restart  analysis  is  specified,  INITAL  reads  the  continuation  data 
and  initializes  the  material  property,  stress,  and  displacement  arrays  for  the 
nodes  and  elements. 

•  SUBSTP.  Subroutine  SUBSTP  controls  the  analysis  of  each  load  case  when 
substeps  are  specified  for  that  load  case.  SUBSTP  divides  the  calculated  equiva¬ 
lent  nodal  point  loads,  applied  displacements,  or  temperature  changes  into  the 
specified  number  of  equal  increments  prior  to  performing  the  analysis. 

•  STRSTF.  The  terms  of  the  global  stiffness  matrix  are  assembled  in  subroutine 
STRSTF  by  sequential  calls  to  QUAD,  JTSTF,  and  BAREL. 

•  OPTSOL.  Subroutine  OPTSOL  solves  the  series  of  simultaneous  equations 
using  Crout  reduction  to  obtain  the  incremental  displacements. 
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STRESS.  Subroutine  STRESS  computes  stresses  and  strains  for  the  two-dimen¬ 
sional  elements  and  prints  the  results.  STRESS  calls  MODCAL,  BAREL,  and 
JSTRES,  used  to  update  the  modulus  values  for  two-dimensional  elements,  inter¬ 
face  stiffnesses,  and  bar  stiffnesses  for  use  in  the  next  iteration  or  load  case. 

QUAD.  Subroutine  QUAD  computes  the  local  element  stiffness  matrix  and 
stress-strain  matrix  for  two-dimensional  elements,  computes  equivalent  nodal 
point  forces  due  to  temperature  changes  of  non-soil  elements,  and  computes 
equivalent  nodal  point  forces  due  to  gravity  forces. 

BAREL.  Subroutine  BAREL  computes  the  stiffness  of  the  bar  elements,  and  up¬ 
dates  the  bar  forces. 

JTSTF.  Subroutine  JTSTF  computes  the  stiffness  of  the  interface  elements. 

JSTRES.  Subroutine  JSTRES  computes  the  stresses  and  relative  displacements 
for  the  interface  elements,  updates  the  interface  stiffness  values,  and  prints  re¬ 
sults. 

MODCAL.  Subroutine  MODCAL  updates  the  modulus  values  assigned  to  the 
soil  elements. 

BUILD.  Subroutine  BUILD  computes  the  nodal  point  loads  which  are  equiva¬ 
lent  to  the  weight  of  the  elements  representing  a  newly  placed  embankment  lift 
and  establishes  the  initial  stresses  and  material  properties  for  these  newly  placed 
elements. 

EXCAV.  Subroutine  EXCAV  calculates  the  stresses  acting  on  an  excavation 
boundary. 

EQNDFO.  Subroutine  EQNDFO  converts  the  stresses  calculated  by  EXCAV  to 
equivalent  nodal  point  forces,  which  are  added  to  the  incremental  load  vector  by 
EXCAV. 

SURFLD.  Subroutine  SURFLD  computes  equivalent  nodal  point  forces  due  to 
a  boundary  pressure  loading  applied  along  the  face  of  an  element  and  adds  these 
computed  forces  to  the  incremental  load  vector. 

SEEP.  Subroutine  SEEP  calculates  equivalent  nodal  point  forces  due  to 
changes  in  the  phreatic  surface  and  adds  these  forces  to  the  incremental  load  vec¬ 
tor.  The  nodal  point  forces  are  formulated  based  upon  changes  in  pore  water 
pressures. 

AUXOUT.  Subroutine  AUXOUT  writes  continuation  data  to  a  file  for  use  in 
subsequent  analyses. 

PRNCIP.  Subroutine  PRNCIP  calculates  principal  stresses  and  the  maximum 
shear  strain  for  two-dimensional  elements. 


•  PRNTFD.  All  non-zero  values  of  the  incremental  load  vector  are  printed  by 
subroutine  PRNTFD. 

•  GETFIL.  Subroutine  GETFIL  initializes  the  execution  of  SOILSTRUCT  by  re¬ 
questing  the  names  of  the  input  and  output  files  and  the  corresponding  opening 
and  closing  of  the  disc  storage  devices. 

•  CNVERT.  Subroutine  CNVERT  creates  a  disc  storage  file  for  use  in  the  plot¬ 
ting  of  computed  results  by  post-processors.  Note  that  CNVERT  discards  all 
air  elements  and  nodes  attached  only  to  air  elements  and  renumbers  the  remain¬ 
ing  existing  nodes  and  elements. 

•  REBAR.  The  compilation  of  the  arrays  for  the  reinforcement  bar  elements  is 
performed  in  subroutine  REBAR. 

•  FILLBARS.  Subroutine  FILLBARS  adds  the  reinforcement  bar  stiffnesses  to 
the  global  stiffness  array. 

•  BARSTIF.  Subroutine  BARSTEF  computes  the  local  reinforcement  bar  stiff¬ 
ness  array. 

•  RESTRESS.  Subroutine  RESTRESS  computes  the  reinforcement  bar  incremen¬ 
tal  strains  and  forces. 

•  NFACTS.  Subroutine  NFACTS  computes  the  shape  function  for  the  reinforce¬ 
ment  bar  elements. 


B3 


APPENDIX  C:  BENDING  OF  STRUCTURAL  MEMBERS 

1.  This  appendix  summarizes  both  simple  beam  theory  and  the  theory  of  elasticity  as 
applied  to  the  bending  of  structural  members.  Simple  beam  theory  is  discussed  in  the 
first  section,  while  the  elastic  response  of  flexural  members  bending  in  plane  stress  and 
plane  strain  are  discussed  in  subsequent  sections. 

Simple  Beam  Theory 

2.  Simple  beam  theory  is  the  most  elementary  of  the  theories  for  flexural  members. 

It  relates  the  flexural  displacements  to  stresses  within  a  bending  beam.  Consider  the 
bending  response  of  the  beam  segment  shown  in  Figure  Cl  (a)  when  subjected  to  a  mo¬ 
ment  of  magnitude  M  (Figure  Cl(c)).  In  a  pure  bending  problem,  the  resultant  force  act¬ 
ing  along  the  x-axis  equals  zero.  Simple  beam  theory  assumes: 

a.  The  material  comprising  the  beam  is  isotropic  and  homogenous. 

b.  Transverse  planes  before  bending  remain  transverse  after  bending;  i.e.,  no 
shear  distortions  take  place. 

c.  The  beam  is  a  straight,  prismatic  member  of  constant  cross  section. 

d.  The  stress  acting  normal  to  the  y-axis  is  zero. 

The  constant  moment  M  causes  the  beam  to  deform  in  the  manner  shown  in  Figure 
Cl(c);  i.e.,  compression  of  the  fibers  along  the  top  of  the  beam  (the  concave  side)  and 
extension  of  the  fibers  along  the  bottom  of  the  beam  (the  convex  side).  The  deforma¬ 
tion  of  the  beam  will  take  the  form  of  a  circular  arc  of  curvature  p  .  Somewhere  be¬ 
tween  the  two  sides,  the  beam  fibers  undergo  no  deformation.  In  this  example,  the 
x-axis  is  located  so  as  to  correspond  to  this  plane  and  is  referred  to  as  the  neutral  axis. 
The  length  of  the  mid-plane  beam  segment  abed  remains  unchanged  after  flexure,  equal 
to  A1  (Figure  C  1(b)  and  Figure  C  1(d)).  Due  to  the  first  assumption,  the  deformation 
varies  linearly  from  the  neutral  axis  to  a  maximum  at  the  extreme  fibers,  as  shown  in 
Figure  Cl(d).  The  elongation  Au  of  the  beam  fibers  at  a  distance  equal  to  -y  below  the 
neutral  axis  is  given  by 

Au  *  -yA0  (Cl) 


with  the  change  in  length  divided  by  the  original  length  equal  to 


Au  -yA8  (C2) 
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Figure  Cl.  Deformation  of  a  beam  segment  according  to  simple  beam  theory 
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To  summarize,  the  derivation  of  the  flexure  formula  using  simple  beam  theory  assumes 
that  the  deformations  within  the  beam  are  either  the  extension  or  compression  of  the 
beam  fibers.  This  theory  makes  no  assumptions  regarding  deformations  in  the  y-axis  di¬ 
rection  nor  in  the  out-of-plane  direction.  This  fact  contrasts  with  the  response  of  a  flex¬ 
ural  member  according  to  elastic  theory,  as  described  in  the  following  two  sections. 

Pure  Bending  of  Beams  in  Plane  Stress 

3.  The  flexural  response  of  beams  bending  in  plane  stress  is  developed  in  this  sec¬ 
tion  using  the  theory  of  elasticity.  The  generalized  Hooke’s  Law  for  an  isotropic,  ho¬ 
mogenous  three-dimensional  solid  is  expressed  in  tensor  notation  as 


1  +  v 
E 


v  c 

E  y  °kk 


(CIO) 


where 

Ejj  -  the  strain  on  the  i  plane  in  the  j  direction 
Gy  =  the  stress  on  the  i  plane  in  the  j  direction 
5jj  =  Kronecker  delta  (  =  1  when  i  =  j  and  =  0  otherwise) 

°kk  =  °1 1  +  °22  +  a33 
Hooke’s  Law  may  also  be  expressed  as 


E 

1  +  v 


+ 


5ij  ekk 


(Cll) 


where  ekk  =  i  +  ^2  +  £33  •  *n  *^e  case  Pure  bending  of  the  Figure  C2  prismatic 
beam  in  plane  stress  (the  out-of-plane  normal  stress  cr33  =  0),  bending  stress  results 
along  the  plane  normal  to  the  Xj-axis.  In  a  plane  stress  problem,  it  is  easily  shown  that 
the  normal  stress  a22  and  the  shear  stresses  o12  ,  o13  ,  and  a23  are  equal  to  zero,  as 
well  as  the  shear  strains  e12,sI3,and  e23  .  The  bending  stress  Ojj  is  expressed  as 


°11  = 


E 

1  +  v 


1  -  v 
l  -  2v  en 


(Cl  2) 
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Equation  C12  shows  that  strains  normal  to  the  x2-axis  and  normal  to  the  x3-axis  contrib¬ 
ute  to  the  value  of  the  bending  stress  Oj  j  in  contrast  to  simple  beam  theory.  From 
Equation  CIO,  the  three  normal  strains  are  given  by 


ell  "  E°H 


(C13) 


*22  =  H°H 


(C14) 


and 


e33  “  E°H 


(Cl  5) 


In  the  plane  stress  case,  there  are  no  restraints  to  movement  in  the  x2  and  x3  direc¬ 
tions,  and  the  strains  normal  to  the  x2-axis  and  the  x3-axis  are  proportional  to  £n  . 


*22  "-^11 


(C16) 


e33  =  -  vel  1 


(C17) 


Introducing  Equations  C16  and  C17  into  C12  results  in  the  stress-strain  relationship 


On  -  E£u 


(C18) 


that  is  equivalent  to  the  relationship  for  the  bending  stress  from  simple  beam  theory 
(Equation  C9). 

4.  The  flexural  response  of  a  bending  member,  modeled  using  a  single  row  of  solid 
finite  elements,  may  be  interpreted  by  using  one-dimensional  bar  elements  as  strain 
gages.  In  this  approach,  as  summarized  in  paragraph  43  for  plane  strain  finite  elements, 
bar  elements  are  placed  along  the  axis  of  bending,  i.e.  parallel  to  the  x  j-axis  ,  and  lo¬ 
cated  at  the  extreme  fibers.  Unlike  simple  beam  theory,  the  e22  and  £33  strains  will 
contribute  to  the  value  of  bending  stress  and  therefore  the  corresponding  stresses  due  to 
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bending  in  a  displacement-based,  plane  stress  finite  element  formulation.  However,  it 
has  been  shown  in  this  section  that  the  value  for  Ojj  in  a  plane  stress  bending  problem 
is  equivalent  to  the  bending  stress  from  simple  beam  theory.  It  will  be  shown  in  the 
next  section  that  this  will  not  be  the  case  when  a  plane  strain  finite  element  model  is 
used. 

Pure  Bending  of  Members  in  Plane  Strain 

5.  The  flexural  response  of  structural  members  bending  in  plane  strain  is  developed 
in  this  section  using  the  theory  of  elasticity.  In  the  case  of  pure  bending  of  the  Figure  C3 
long  member  in  plane  strain  (the  out-of-plane  normal  strain  E33  =  0),  bending  stresses 
result  along  the  plane  normal  to  the  Xj-axis.  In  this  problem,  the  normal  stress  o22  is 
equal  to  zero,  as  well  as  all  shear  stresses  and  shear  strains.  O33  is  not  equal  to  zero,  un¬ 
like  the  plane  stress  problem.  The  bending  stress  Oj  j  is  expressed  as 

°n  =  (TTvjo  -  2v)[(1  “  V)£1I  +  ve22]  (C19> 

The  strain  e22  can  be  removed  from  Equation  Cl 9  because,  as  for  the  plane  stress 
case,  there  is  no  restraint  to  movement  in  the  x2  direction.  To  accomplish  this,  the 
stress  boundary  conditions  (o22  =  o33  *  0)  are  introduced  to  Equation  CIO.  The  result¬ 
ing  normal  strains  are  given  by 

o  _  (1  +  v)(l  -  v)_  (C20) 

ell  E  n 


and 


(1  +  v)  v  , 

E  °H 


(C21) 


Strain  normal  to  the  x2-axis  ,  expressed  in  terms  of  the  strain  normal  to  the  Xj-axis  ,  is 
given  by 


*22  ~  “ 


v 

1  -  veH 


(C22) 
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By  combining  terms.  Equation  C20  becomes 


(1-v2)  11 


(C23) 


By  following  the  same  steps  described  on  page  C3,  the  moment  M  is  expressed  as 


M-— SL-Si 
(i  -  v2)  y 


(C24) 


where  y  =  distance  from  the  neutral  axis  to  the  beam  fiber  at  which  the  pure  bending 
strain  etl  is  measured.  Equation  C23  differs  from  the  relationship  for  bending  stress 


*3 

(OUT-OF -PLANE) 


Figure  C3.  Pure  bending  of  a  long  member  inplane  strain  *  0) 
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in  simple  beam  theory  (Equation  C9)  by  the  factor  1/(1  -  v2).  When  a  single  row  of 
solid  finite  elements  are  used  to  model  the  member  bending  in  plane  strain,  one-dimen¬ 
sional  bar  elements  are  located  along  the  axis  of  bending,  i.e.,  parallel  to  the  Xj-axis. 
The  622  strain  will  contribute  to  the  value  of  bending  stress  and  therefore  the  corre¬ 
sponding  moment,  as  discussed  in  paragraph  42.  Thus,  when  computing  the  internal 
moment  that  is  equivalent  to  the  a  j  j  distribution  across  the  section,  using  the  one-di¬ 
mensional  bar  strains  £j  j  from  the  plane  strain  finite  element  analysis,  a  factor  equal 
to  1/(1  -  v2)  is  included  in  the  relationship  (see  Equation  31). 


APPENDIX  D:  NOTATION 


be  Unit  width  of  sheet-pile  wall  mesh 
B  Bulk  modulus 
c  Cohesion  intercept 

Cj  Cohesion  intercept  for  interface  elements 
d  Unit  thickness  of  slice 
E  Young’s  modulus 
Ee  Equivalent  Young’s  modulus 
Ej  Initial  Young’s  modulus 
E(  Tangent  Young’s  modulus 
Eur  Young’s  modulus  during  unloading  or  reloading 
Ej  Vertical  stiffness 
E2  Hoop  stiffness 

[F]  Incremental  load  vector  applied  to  the  nodal  points 
G2  Shear  modulus 
H  Vertical  distance 
I  Moment  of  inertia 
kn  Interface  normal  stiffness 
ks  Interface  shear  stiffness 
ksi  Initial  interface  shear  stiffness 
k$t  Tangent  interface  shear  stiffness 
K  Modulus  number 

K  Constant  relating  undrained  Young’s  modulus  to  undrained 
strength 

K  Uniformly  distributed  radial  spring  acting  to  restrain  radial  move 
ment 

Kb  Bulk  modulus 

[Kb]  Stiffness  for  reinforcement  bar  element 
Kj  Interface  modulus  number 
K0  At-rest  earth  pressure  coefficient 
Ks  Planar  spring  stiffness 
Kur  Unload-reload  modulus  number 
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[K]  Global  stiffness  matrix 

m  Bulk  modulus  exponent 

M  Moment  per  unit  width  of  wall 

n  Modulus  exponent 

n;  Interface  modulus  exponent 

p  Lateral  pressure  inside  cell 

p'c  Effective  consolidation  pressure 

P  Radial  pressure  inside  cell 

Pa  Atmospheric  pressure 

rc  Radius  of  curvature  for  bending 

R  Cell  radius 

Rf  Failure  ratio 

^  | 

Rfl  Failure  ratio  for  interface  element 

I 

Su  Undrained  strength 
SL  Stress  level 

SL|  Stress  level  for  interface  element 

T  Unit  interlock  force 

i 

T.  Spiring  force 
s  r 

t  Cell-wall  thickness 
t  Slice  thickness  (=1) 

u  Node  displacements,  radial  displacement  of  nodes 
u  Incremental  nodal  displacements 
ub  Displacements  of  bar  element  nodes 
A,,  Nodal  relative  interface  elemnt  displacements 
As  Average  relative  interface  element  shear  displacements 
AOj ,  A C2 .  Ao3  Changes  in  values  of  principal  stress 

5  Angle  of  interfal  friction  along  the  interface 
S  lateral  movement  along  the  wall 
ea  Axial  strain 
Ej,  Bending  strain 
Zy  Change  in  volumnetric  strain 


e, ,  £j  Outer  fiber  strains 
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V  Poisson’s  ratio 
Vj  Unfailed  Poisson’s  ratio 
Yw  Unit  weight  of  water 
Gj  Major  principal  stress 
g'j  Major  principal  effective  stress 
G3  Minor  principal  stress 
g'3  Minor  principal  effective  stress 
G  jj  Horizontal  effective  stress 
Gy  Vertical  effective  stress 
Gn  Normal  stress 
ot  Tensile  strength 
(Gj  -  g3)  Current  deviator  stress 
(Gj  -  G3)f  Deviator  stress  at  failure 
(Gj  -  G3)ult  Ultimate  deviator  stress 

t  Shear  stress  -  interface  element 
Tf  Shear  stress  at  failure  -  interface  element 
xult  Ultimate  shear  stress  -  interface  element 
9  Angle  of  internal  friction 
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